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SUMMARY 


Forced  response  calculations  are  performed  on  the  H-21  helicopter  equipped 
with  floating  fuel  wings  for  the  0%  and  100%  fuel  configurations.  Using 
third  harmonic  shaft  loads  measured  on  an  H-21  helicopter,  and  a  deflection 
test  obtained  analytical  model  of  the  H-21  fuselage,  the  Associated  Matrix 
Method  is  employed  to  calculate  the  fuselage  forced  response.  Solutions  are 
obtained  on  the  Remington  Rand  1103A  computer  installation  at  Wright  Field. 
Results  of  the  response  calculations  are  presented  as  modal  time  histories 
of  the  fuselage  for  third  harmonic  azimuth  positions  of  0,  45,  90,  and  135 
degrees,  and  the  resultant  vertical  and  lateral  responses  at  the  cockpit 
floor  during  an  rpm  sweep. 

Results  of  the  analysis  indicate  that  at  the  0%  fuel  configuration,  vertical 
response  compares  closely  to  that  of  the  standard  helicopter.  Laterally,  the 
response  increases  with  rotor  speed,  having  an  acceptable  level  at  the  lower 
rotor  speed  range  and  somewhat  higher  amplitude  at  268  rpm.  For  the  100% 
fuel  configuration,  the  forced  response  calculations  indicate  the  presence  of 
a  lateral  natural  frequency  close  to,  but  below,  the  3_f2-  excitation  at  250  rpm 
In  both  the  lateral  and  vertical  directions,  the  calculated  undamped  vibration 
levels  are  high  at  the  lower  rotor  speeds,  but  appear  satisfactory  at  268  rpm. 
Past  experience  indicates  that  predicted  vibration  levels  near  resonance  are 
reduced  by  damping.  Nevertheless,  as  a  result  of  the  calculations,  it  is 
recommended  that  a  simplified  ground  shake  test  be  performed  to  substantiate 
the  calculations  and  provide  some  reasonable  prediction  of  the  effect  of 
damping. 
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SECTION  I 


INTRODUCTION 


The  present  wing-fuselage  forced  vibration  study  is  part  of  an  analytical  and 
wind  tunnel  study  being  conducted  under  the  Reference  1  Transportation  Research 
Command  contract.  The  program  is  aimed  at  the  development  of  a  means  for 
helicopter  ferry  range  extension  through  application  of  a  floating  fuel  wing 
concept.  Ap  initial  feasibility  study  of  the  floating  wing  concept  was  con¬ 
ducted  by  Vertol  under  an  earlier  contract,  Reference  2,  and  the  results 
reported  in  Reference  3.  The  present  analytical  and  wind  tunnel  investigation 
is  under  the  Reference  1  contract  and  is  based  on  a  Vertol  proposal,  Reference 
4.  This  yibration  investigation  is  a  part  of  Phase  II  of  the  contract.  Addi¬ 
tional  dynamic  studies  of  the  contract  include  ground  and  air  instability  re¬ 
ported  in  Reference  5  and  flutter  calculations  reported  in  Reference  6. 

The  range  of  present  helicopters  with  normal  fuel  load  is  less  than  400  nau¬ 
tical  miles.  Even  with  additional  internal  tanks,  the  helicopter  range  is 
less  than  11Q0  miles.  With  floating  wings,  the  range  can  be  extended  to  as 
much  as  2400  miles,  corresponding  to  the  longest  over-water  distance  on  the 
Pacific  Ocean  ferry  route. 

Each  floating  wing  contains  compartmented  fuel  tankage  connected  by  lines  to 
the  helicopter's  main  tank.  The  wing  lift  supports  the  fuel  weight  that  it 
carries,  and  the  helicopter  acts  as  a  tow  to  propel  the  wing  forward.  Wing 
attachment  to  the  helicopter  is  through  a  hinge  so  as  to  eliminate  the  bend¬ 
ing  moments  applied  at  the  fuselage  by  conventional  wings,  thus  avoiding  the 
addition  of  extensive  wing  carry  through  structure  to  the  helicopter.  The 
hinge  line  is  not  longitudinal,  but  is  skewed  aft  as  shown. 
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As  fuel  is  consumed  and  the  wing  becomes  lighter,  it  tends  to  flap  upward 
about  its  hinge.  Because  of  the  skewed  orientation,  the  angle  of  attack 
at  any  chord  line  is  reduced  as  the  wing  flaps  up,  the  lift  is  reduced,  and 
the  wing  assumes  a  new  mean  position.  Full  span  pilot  controlled  wing 
flaps  are  also  provided,  so  that  the  trim  attitude  of  the  wing  may  be 

adjusted;  these  are  also  used  as  high  lift  devices  during  the  running 

takeoff . 

In-flight,  the  wing-helicopter  combination  is  excited  by  rotor  induced 
loads  such  as  measured  on  the  H-21  helicopter  under  Reference  7.  In  a 
previous  study  performed  under  Reference  8,  the  forced  response  of  the 
H-21  helicopter  without  floating  fuel  wings  was  determined  theoretically 
by  applying  measured  shaft  loads  to  an  analytical  representation  of  the 
fuselage  which  included  vertical-lai_eral  coupling.  The  present  analyses 
use  the  coupled  fuselage  program  of  Reference  8  programmed  on  the  Remington 

Rand  1103A  computer.  To  include  the  floating  fuel  wings,  the  H-21  fuselage 

representation  is  expanded  to  include  an  effective  wing  matrix  which  repre¬ 
sents  the  flap-torsion  wing  modes. 

The  fuselage-wing  forced  response  is  investigated  here  for  the  wing  empty  - 
0 %  fuel  and  wing  full  -  100%  fuel  at  each  of  the  possible  operating  rotor 
speeds  of  250,258  and  268  RPM  for  a  cruise  speed  of  90  knots. 

Section  II  presents  the  method  of  analysis,  followed  by  a  discussion  of  the 
results  in  Section  III  including  a  description  of  the  analytical  represen¬ 
tation.  Conclusions  are  given  in  Section  IV,  and  References  in  Section  V. 
The  matrix  derivations  are  presented  in  Appendix  A,  along  with  the  method 
of  solution  in  Appendix  B. 


SECTION  II 


METHOD  OF  ANALYSIS 


A.  General 


The  method  used  here  for  forced  mode  calculations  is  the  coupled  vertical- 
lateral  analysis  of  Reference  8.  A  coupled  vertical-lateral  system  was 
chosen  for  representing  the  H-21  fuselage  as  a  result  of  previous  investi¬ 
gations,  both  analytical  and  test  indicating  the  significance  of  coupling. 

The  structural  deflection  test  of  the  H-21  fuselage  under  Reference  9,  and 
the  initially  limited  correlation  shown  between  uncoupled  analytical  and 
ground  shake  test  results  in  Reference  10  indicated  that  an  analytical  model 
of  the  H-21  fuselage  must  include  vertical-lateral  coupling. 

The  analysis  includes  the  following  basic  matrix  types:  (1)  discrete  lumped 
masses  and  moments  of  inertia,  (2)  suspended  lumped  masses  and  mass  moments 
of  inertia,  (3)  bends  permitting  the  selected  elastic  beam  axis  to  follow  the 
general  fuselage  shape,  (4)  sections  of  weightless  elastic  beam  which  permit 
the  inclusion  of  bending,  shear,  tension  and  torsional  stiffness  properties, 

(5)  concentrated  springs  which  simulate  the  attachment  of  large  local  mass 
and  inertia  items  to  fuselage  structure,  such  as  the  rotor  transmission,  and 

(6)  ground  spring  matrices  which  permit  attachment  of  the  vibrating  body  to 
an  external  ground  as  in  a  bungee  suspension  for  ground  shake  testing.  Also, 
to  provide  the  specialized  representation  necessary  for  the  H-21  fuselage  and 
to  include  forced  response  calculations,  the  above  type  of  matrices  were  ex¬ 
panded  to  include  (1)  sections  of  weightless  elastic  beam,  which  permit  frame 
distortion  in  combination  with  bending,  shear  and  torsion  stiffness  properties, 
(2)  stations  which  permit  the  attachment  of  a  coupled  vibration  system  such  as 
the  flexibly  mounted  engine  with  known  modes  and  frequencies  and  (3)  force 
matrices  which  permit  exciting  forces  to  be  applied  at  any  point  on  the  vi¬ 
brating  body,  such  as  the  external  shaker  in  a  ground  shake  test  or  the 
measured  flight  loads  at  the  rotor. 

This  existing  program  is  adapted  to  the  H-21  fuselage-wing  combination  using 
an  effective  wing  matrix  to  transmit  the  dynamic  wing  loads  to  the  fuselage 
with  the  uncoupled  wing  modes  as  generalized  coordinates.  As  derived  previ¬ 
ously,  with  the  addition  of  the  effective  wing  matrix,  the  matrices  and  ele¬ 
ments  are  presented  in  Appendix  A  along  with  the  method  of  solution  in 
Appendix  B.  Figure  1  presents  the  digital  computer  flow  chart  showing  the 
basic  programming  procedure  for  the  Remington  Rand  Univac  1103A  digital  com¬ 
puter.  Application  to  the  wing-fuselage  matrix  program  is  illustrated  by 
Figure  2,  Computer  Input  Forms,  partially  completed  for  the  H-21  helicopter 
equipped  with  floating  fuel  wings.  Figure  3,  Input  Sequence  Chart  for  the 
Coupled  Matrix  Program  and  Figure  4,  Numerical  Output  for  a  Fuselage-Wing 
Forced  Mode. 
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FIGURE  1 .  COUPLED  VERTICAL- LATERAL  MATRIX  FLOW  DIAGRAM 
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Rotor  Shaft  Loads 


INPUT  FOR  H~21  HELICOPTER  EQUIPPED  WITH  FLOATING  FUEL  WING 


RQ,  X  051 

HI, 

032861, 

1021 

27100, 

352.00, 

0, 

90, 

250, 

1, 

R35SF  2,50  x  102 
X  1, 

50, 

1, 

B  10, 

X  000 
78.6, 


Case 

Job 

Date 

Aircraft  Model  No. 

Gross  Weight 
C.  G.  Station 
Flight  No. 

Airspeed  Knots 
Rotor  RPM 

Engine  (0  -  Rigid,  1  -  Suspended) 

Blades  (l-On,  O-Removed) 

Rotor  Speed  RPM 

Collapsed  Matrix  (0  -  No  Print  Out,  1  -  Print  Out) 
Total  Matrices 

Program  (0  -  Free,  1  -  Forced) 

Boundaries  (10  -  Free  -  Free) 

Number  of  Forcing  Frequencies 
Forcing  Frequency 


MATRIX  INPUT 
(Read  Down  Columns) 


FIGURE  2  COMPUTER  INPUT  FORM 


6 


Transmission  Transmission  Elastic  Bay  Rotor  Shaft 

Mass  _  Spring 


Loc , 


Code  &  Elements 


+1.05  x .  10 ' 
+4.00  x  106 


B07006 


B01015 

0 

0 

0 

0 

1.56  x 
1.00  x 
1.0 
1.0 
1.0 
1.0 
1.0 


20 


10J 

10 

1020 

1020 

1020 

1020 

1020 


E 

G 


F+l .  0 

X 

1020 

W« 

W 

C 

+1.0 

X 

1020 

Kz 

td 

+1.0 

X 

1020 

6 

Kck 

H 

+2.5 

X 

107 

K/* 

+2.5 

X 

107 

Kjf 

L 


■ 

+1.05  x  10 ‘ 
+4.00  x  106 

E 

G 

5 

B05007 

1.0549  x  10° 

M 

5.6  x  101 

u 

6.1  x  101 

6.1  x  101 

I* 

-3.0  x  10° 

a 

0 

b 

0 

c 

c 

o 


vi  bO 

'  C 


u 

CL 


<d 

PQ 


tn 

cd 


W 


d 

o 


CL, 


d  xj 
o  d 

T—t  QJ 

>>  CQ 
P+ 


Loc  . 

Code  & 

Elements 

6 

B07006 

F+l  .  0  x 

1020 

Kx 

+3.5  x 

105 

Ky 

+2.3  x 

104 

*<z 

+1.0  x 

1020 

+4.0  x 

107 

+7  0  x 

107 

7 

BO 10 15 

F  0 

b 

0 

"e 

0 

tan 

0 

cn 

+3.05  x 

10 1 

L 

1 .0  x 

1020 

Ax 

1.0  x 

1020 

Ay 

1.0  x 

1020 

A> 

1  0  x 

1020 

Jx 

1.0  x 

1020 

Iv 

1.0  x 

1010 

-1 

*  z 

105  x 

10' 

E 

4.00  x 

106 

G 

8 

B03003 

F+l. 6505  x  10" | 

Cos  P 

-9.8629  x  lO’1 

Sin/S 

+2.795 

x  10^ 

P 

FIGURE  2  (CONTINUED)  COMPUTER  INPUT  FORM 


TYPE 

NPUT  CODE  -  INPUT  SEQUENCE 


1 


INPUT  SEQUENCE  FOR  COUPLED  MATRIX  PROGRAM 


Lu 

O 


r>:  m 

y  o 
2  o 
<£ 


lu 

o 


E-» 

to 


CL 

K  o 

r  in 

O  nj 


s 

a 


Ld  o 


O 
CC  * 
Id  co 


iJ  O 


.0 


lu  m; 

►-  to 
<C  «M 
O  9 


Pi 

o 

Cjh 


o 

IT  3 

ip  o 

CL  O 

o 


•  o 


1  CO 


O  d(J)  ^  CM  CM  AJ  f"  —4  4 

>  4  4  >  J  j  p  r\J  o  J 

‘  ■  -4  *n  4  4  «n  >  >n 


I  — •  r-  r-  O'  o  O'  O'  ^  to  r-  <M  o  cn  us  c-~  ns  :foOwo'g'rj^e>iM  '-5 

- -  -  -  -  -~  -  i-  -  |  q  _y  -f  <S  In  o  Is  O  .5 


j'oooo  n  im  >  ft  n  p  -t  «  >  in  N  Li  in 


n  (0  a-  o  o  m  r  -  t*' 

3  P  a|  OOO  -*  -* 


VJ 

£ 

O 

•o 

M  AJ  -4  4  t  CD 

t  1  1  1  1  1 

4  n  *1  n  i/t 

1111 

4  -• 

-a  _i  — «  — <  0  m  i-t  »n  0  —* 

1  1  1  1  1  1  1  1  1 

•  •  *  • 

- 1  —4  •— 4  —4 

IIII 

1  1 

•  • 

1  1 

•  • 

n  4 

1  i 

•  •  4  • 

n  s4  in  aj 

1111 

«  • 

VM  P 

1  1 

P 

1 

1 

P 

1 

iD 

1 

P 

1 

4  * 

1  1 

CM 

1 

- 

T 

OJ 

p  ^  r-1  ^  r^i  -f 

t-/4  -rt 

4  vit  »n  Lit  t  4  -4  -4  r’t 

rati-t  t'4UtininLtutLtLtutJo*n,mioioin  104-4 

-4 

4 

<T 

J4Tin 

J  O  O  O  O  )  o  o  o  o  OO  O  O  O  o  O  O  O  3  3 

O  30  0  0  3  3  3 

3  3 

3  D 

30 

300000 

O 

O 

O 

O 

030 

O 

y 

& 

o 

n>  r- 

4  hh,s 

in  rc  cn 

4  — 

x  id  3  -  0  n  0 

C\  O'  >  — 

ut  r- 

— <  ft 

ft  d 

4 

ft  >  -0 

>  AJ 

P 

P 

X> 

P 

p 

us  in 

O 

4 

J 

«M  t»  00 

hn^  in 

in  4  Ln  fn  cS 

*t  f- 

in  V-  co  4  0  ft  -< 

C3  O  a)  co  d-  O'  co  ut  im  -4 

t - - 

r-  ut  -4  10 

—4  P 

Tf 

-f  ^t  £3 

& 

CT' 

\n 

*-  r» 

o  >  O'  in 

In  >  4  ut 

<M  Q 

3  n  Aimco<noco 

0  0  U3  UD 

OQh  H 

■  4  O' 

O' 

c\J  M  4 

4  SO  P 

p  'n 

n 

3 

4  -r 

n  >it 

i  4  ' 

t  i 

till 

11  1  1  1  r  j 

1  »  1  1  1  1  i  t  1 

‘‘ 

r 

n 

rt  *n 

ft  -n  ft  m 

*n  4  t  -1  1 

■n  -4  4  •n  cn  ft  m 

4  4  4  4 

4  4 

4  -4 

4 

4  4  4 

4  4 

4 

ft 

•n 

ft 

m 

rn  m 

ft 

— 

c 

c  c 

C  O  O  o 

0  0  0  c?  c: 

C  C 

c  0  0  0  0  0  0  0  C  C  C 

c  c  0  0 

0  c 

c  0 

c  0 

r~ 

0  c  0 

C  C 

Li 

c: 

C' 

0 

T 

£0 

in  ut  m  Via  ^  ^  n 

>00  i.n  in  \n  ut 

-•fro  A  r-  Al  M  0  L3  -o  03 

0  O'  0  <0 

O  -4 

—4  Ut 

-O' 

<M 

OOO  AP'O'^SJ  'O 

-4 

—4  -1 

-4 

-4 

-J 

in 

in  in 

in  'O'O'O 

C\J  « 

m  m 

in  M3  so  co  <0  co  (0 

vO  '0  '0  AJ 

AJ  4  4  r- 

r-  a? 

CO  M3 

r- 

4 

POO 

O'  CT' 

—4  —4 

4 

>n  in  w  i 

tn  n  t/i  in 

g3  ^  ^  _H 

OOSONCM^  4  4 

rt  -4  — »  — ♦ 

"4  'H  -4  —* 

~t  -4 

-4  -4 

4  m  0 

-4  4-4 

O  O 

0  at  00  <30  <T>  O' 

0 

in 

* 

o 

ft 

ft  n 

•n  Ann 

1  1  1  1  1 

1  1 

4  13  -+  co  MO  -4  -< 
II  If 

CM  —•  —•  CM 

1  1 

CM  OJ  CM  Al 

CM  <M 

M  M  i\J  ^ 

—4 

_ 

—4  ^ 

- 

tn  n  ft  rt  h  n  *n 

1  t  1  1  1  1  1 

- 

u 

M  \J 

M  HIM  M 

AJ  CM  TM  r>l  AJ 

n  n  net 

ft  ft 

n  ^i  fn  n  n  it 

m  m 

m  fn 

n 

n 

\! 

M  M 

ai  m 

3 

3  O 

O  J  J  J 

OO  3  0  0 

-3  '  3 

O  .3  3  .3  .3  0  0 

0  0  3  3 

3  3  0  0 

0  0 

0  0 

0  3 

3  .O  O  O 

3  3 

00 

3 

3 

O 

3  O 

0 

V 

-J 

3 

3  3 

olnoinlno'O'44 

m  O  i>  -• 

©nncti 

N  h-  4 

O 

0  in 

U  M 

p 

-4 

-4 

—4 

O 

O  O 

9 

-4 

— * 

—4  -  < 

<  — «  — t  *-+ 

co  r-  '4'jvjrvigin--'-'^A,'' 

CM  -H  OO  <D  O  O 

~4\p 

Q 

0  ut  Ut 

Lt  lit 

t 

4 

t 

4 

LS  CO  Gki 

CO 

O' 

N 

r- 

r-  i'- 

0  n  0  -4 

■  -* 

O  r>  -t  ft  ft  -4-4 

4  A/  AJ  if) 

'0  r-  r-  co 

co  O' 

.>  -> 

n  in 

O 

0  C3  .  0 

49  p 

r- 

3 

r3 

O 

P 

co  <n 

•-0 

O 

o 

aj 

i 

CM  CM 

1  1 

iM  4t  n  n 

iiii 

!M  J)  O  D-0 

t  1  1  1  1 

4  4 

1  1 

-•  -4  ft  ft  ft  -4-4 

1  1  1  1  1 

-a  -4  -4  AJ 

AJ  Al  AJ  AJ  <M  Al 

AJ  AJ 

n  -« 

1 

AJ 

1 

rg  -4  -4 

1  1  1 

• — 4  »4 

1  1 

p 

—4 

1 

1 

P 

1 

■0  D 

1  1 

p 

l 

—4 

Al 

AJ  CM 

AI  AJ  Al  AJ 

AJ  -*  AJ  AJ 

AJ  M 

-4-4-4 

-4 

-4  -* 

AJ 

AJ  M  M 

N  Ai 

Al 

M 

Aj 

AJ 

AJ 

M  AJ 

AJ 

■n 

£C 

O 

O  O 

o  o  o  o 

OOO  30 

3  3 

O  0  0  3  0  0  0 

0  0-30 

3  O  -3 

o 

O  O 

O  3 

O 

D  3  O 

o  0 

O 

O 

O 

3 

3 

3  O 

0 

V 

0 

3  5 

t  -•  -< 

^  AJ  J 

frnwiKI  ;J  ill  ill  h.M 

OO  3  > 

0-4-4  0 

0  03 

O  O 

0  > 

O' 

->  0  0 

O  3 

O 

3 

3 

3 

3 

3  O 

3 

> 

4 

4  4 

4  4  4  4 

4  -\J  M  ft  ft 

•0  0 

M  M  vO  >0  >0  '0  0 

ut  0  m  3 

3  3  i"N 

n-  it 

it  it 

il  -* 

-4 

— 4  r~  S 

—4  —4 

O' 

O' 

4 

.M 

Al  AJ 

AJ 

f- 

>- 

4 

4  4 

4  ft  ft  n 

ft  .-0  03  -4  -* 

4)  41 

44  -4  —4  —4  at  at 

0  0  0  4 

4  a-  r~  0 

>  -4 

-4  4 

4 

4  ft  n 

n  n 

Ai 

AJ 

O 

O  O 

O 

f- 

O 

- 

-  - 

-♦  •  ■<  H  •< 

-f  V/>  M3  -* 

IIII 

1  1 

MMlOOiOM^NNK  4 

1  1  1  1  1  I  1  1  1  1  1 

4  -^4  4  in 

1  1  t 

IO  -4 

-  - 

-*rnrnrnmmtncn 

1  1  >  1  1  1  1 

AJ  CM  AJ  CM 

1  1  1  r 

1 

1  1 

1 

1 

_ 

-4 

rg  ,M  <M 

Al  CM  rj  vj  <M  AJ  i~J 

(\J  AJ  -4  “4  -4  VJ  'J 

CM  CM  AJ  cn 

»n  t-t  ri  Kt 

ft  ct 

m  m 

cm  «n 

CJ  AJ  CM 

N  CM 

AJ 

A!  c\l  AJ 

At 

O 

O  O 

A  0  5  O 

'3  0  0  0  O 

O  O 

3  3  0  0  O  O  O 

0000 

3  0  0  0 

3  0 

0  3 

0  0 

O 

3  3  0 

0  0 

3 

O 

'3 

y 

_J 

a 

o  o 

Q  AJ  AJ  AJ 

a-  n  n  O'  O' 

ift  iTt 

o  i/t  it  M3 

M3  iO  iO  r>* 

r--  —i 

-4 

r-  in 

CO 

X)  4  4 

P  SO 

n 

-4 

-4 

-4 

O 

O  O 

O 

0 

p 

O  o 

a  o  o  o 

o-«-«oo 

Al  J 

•n  n  -•  -4  -•  -•  4 

4  0  0>M 

AJ  a)  t)  n 

n  0 

O  0 

•3 

OOO 

0  0  CM 

it 

n  0 

3  3 

•> 

CO 

X 

> 

O'  > 

X  n  n  n 

^  CO  <3  N  J--  PJ  M 

>>101010  00 

>0.  i.-i 

‘OlftdtM3'Ovaso\QO  4 

-4  a-  r-  r-  n- 

D 

P 

p 

O 

9  0 

3  «5 

o 

-o 

O  0 

o  -H  -1 

4  «t  fi  nnnn 

1  1  1  1  1  !  1 

rJ  aj  \0  VO  m3  -4  -4 

-•  O  X)  -• 

-4  ,4  -4 

.4  —4 

-  -1 

it  -* 

cX3 

POO 

CD  CD 

0 

•0 

l 

P 

1 

p 

1 

0 

O  '•3 

•O 

*M 

o 

4 

4 

4  t 

4  4  4  4 

■n  m  n  n  rt 

Jt  O 

•  O  .it  at  O  i  t  .it  at 

•.n  in  0  n 

0  in  in  -o 

in  10 

iO  n 

tn  it 

•n 

n  m  n 

n  .-4 

4 

n 

n 

4 

4 

4  4 

4 

4 

§ 

O 

O 

O  O 

9  9  9  9 

0  0  c-  0  v 

0  9 

C  O  C  C  Y  c>  ^ 

1'  C  O  13 

9  9  9  p 

p  0 

9  9 

<3  0 

0 

9  9  9 

0  0 

0 

0 

0 

O 

0-  CO 

O 

O 

4 

4 

3 

i  .-i  -!  fi. 

n-  r-  ii  at  0 

6  4 

4  4  4  j-  -O  d  »n 

n  .n  ai  iM 

O  O  AJ  (VJ 

fT- t 

n 

n  in  .\i 

.M  > 

O' 

.0 

ip 

-<  i> 

.1 

.-4 

<r 

p 

P 

O  4 

4  J  oo 

0  O  -J  3  «c 

O  sC 

c  'C  'C  a:  m;  'O  -n 

•lO  if  A  h 

in  in  —  — 

O  M3 

A-  A- 

in.  in 

r-  r-  r- 

f--  P 

p 

A- 

r-~ 

p 

P 

P  M 

ft 

ft 

1 

< 

U 

(>i 

Aj 

O-  U 

O  O  4  4 

0  id  f*n  >o 

O  -Ni 

l\  l'<  l\|  ill  O'  l>  -4 

-•  0  h-  A- 

0  0  dy  > 

ft  ft 

O  u 

10  CO 

.3 

0  O  O' 

4  P 

p 

-0 

p 

p 

MJ 

P  4 

-* 

y 

o 

1 

•n 

t 

‘M  M 

1  I 

M  a  — ♦  — < 

iiii 

j>  at  n  at  Jt 

1  1  1  1  1 

at  0 

1  1 

0  0)  M3  M3  =r  4  <> 

IIII 

O'  m3  4  4 

1  1  1 

At  r*t  <AJ  M 

IIII 

>M  IM 

1  1 

.AJ  AJ 

1  1 

so  MJ 

M 

1 

M  M  4 

1  1 

4  -« 

> 

ft 

n 

n  lit 

A- 

9 

4 

4 

4  4 

4  4  4  4 

4  44-  4  4 

J-  4-444-4UtCtiO 

Vt  to  lit  4-0 

3  4  4  *4  4-  4 

-r-r- 

4- -4 

-4  4  "4  -r  -4  -4-  >4 

4  <f 

n 

•t 

n  n 

m 

n 

O 

O 

3 

3  3 

3QOO 

3  0  0  0  0 

O  3 

O  .3  3  3  O  3  O 

0  3  O  O 

OOO  O 

3  0 

0  0 

3  O 

O 

O  3  O 

O  3 

p 

O 

0 

3 

0 

0  3 

0 

0 

n 

n 

O'  4 

4  4^-^ 

-4-4  0-0  Al 

(M  -O 

i/t  at  -O  lO  1 M  M> 

•0  1 6  ,n  n 

JO  lO 

'  0  0 

m3  '6 

'0 

M3  '6  6 

p  r- 

A- 

A* 

0 

0 

6  CQ 

p 

p 

A  t 

J  f  rt  n 

ft  n  0  0  \J 

J  J 

''J  \J  J  M  4  f  ^ 

'OOt'*' 

ut  0 

it  in 

in  ut  in  > 

O  co 

O 

0  cq  in 

Jt 

Ut  4  in 

Lt 

< 

m 

V 

-•  M 

j  m  a.  n* 

>  4  4  0 

0  * 

^  >  >  >  4  t  ^ 

'  a  i<t  n 

0  0  in  in  q3  co 

-4  -4 

.  -4 

•* 

-4  -4  CM 

M  P 

0 

0  M3 

0 

p 

p  rsi 

p 

P 

.d 

CD 

O' 

i 

O' 

i 

0  P 

1  1 

•o  o  -t  -n 

iiii 

n  n  t  n  ai 

1  1  *  1  1 

CM  -• 

f  1 

— 4—4  '4  — <fN.  A»4 

4  j-  0  0 

1111 

0  ft 

it  m 

l  1 

.0  in 

1  1 

nim/i  .0 

1111 

P  fM. 

1  1 

A- 

1 

r- 

1 

'  1 

1 

1 

-4  AJ 

1  1 

AJ 

1 

AJ 

I 

a 

4 

4 

in  n 

>n  ifi  ut  n 

4  4  4  4  0 

0  n 

n  'O  ut  n  in  04 

4  r  ./-  4 

4  4 

4  4 

4 

4  4  4 

t  P> 

in- 

t 

4 

4 

4 

4  -4 

4 

4 

< 

O 

O 

O  o 

9000 

O  O  O  O  O 

O  0 

000c  0  0  0 

O  O  <3  O 

0000 

0  0 

p  p 

O  O 

0 

p  9  9 

9  ? 

C3 

O 

0 

0 

y 

O  O 

O 

C 

DC 

O' 

> 

4  4 

4444 

>  l>  O  O  -4 

J-,  i, 

--4  '4  -4  0  it  rj 

O  0  '-A 

Al  M 

A-  A- 

A- 

A-  it 

.n  n 

n 

Jt 

i/t 

at 

n 

It  it 

P 

O 

a 

o 

O 

■n  ft 

•*t  n  *n  n 

■  M  rj  -4  -<  n 

n  -j 

j-  4  4  4  -4  -4  4 

4  GO  'J  O 

dUHJOO'it 

>  4 

3)  O 

> 

4>  O'  4 

.t  —4 

-4 

O' 

3 

> 

> 

4 

> 

CV 

(O 

at  at 

U1  dt  ifl  Ul 

ft  ft  ft  ft  --4 

'4  O 

O'  0  i>  4  lO  'O  — 4 

-4  —  n  ;n 

"4  —4  — 4  — ' 

O'  O' 

ft  m 

ft  ft 

ft 

ft  ft  4 

3-  CM 

v-\J 

J 

4 

j 

4 

4  4 

-* 

-J 

< 

AJ 

<AJ 

'A  ft 

ft  /n  <n  n 

-4  — «  .-i  ,H  AJ 

IIII 

rg  -4 

-•  -4  -•  -4  rj  CM  4 

4  -r  4  4 

4  4  4  4 

AJ  .M 

Al  iVl 

AJ  OJ 

'M 

<AI  JSJ  -* 

-4  P 

1 

p 

1 

- 

-4 

-1 

-4 

.-4  -4 

4 

4 

Al 

AJ 

M  v 

M  M  A  >n 

•a  m  n  -nil 

M 

aj  M  M 

Al  -4 

Al 

AJ 

;ai 

M 

M  AJ 

3 

- 

3  ■? 

'  f  (- 

3  3  3  3) 

:  0 

'  ;  0  :  i  ;  3 

3  '  3  3 

3  3 

3  3 

3  3 

3 

' 

3 

,3  j  ’ 

3 

1. 

r-  r1- 

ft  /'l  on  r~ 

O  O  O  O  Ul 

...  .> 

i  i:  d  i 

fwn 

AJ  O' 

O' 

iM  Aj  —4 

-4  r  - 

A- 

P 

(M 

^-4 

-4 

A-  [9 

vM 

Al 

4 

4 

p—  n 

<t  4-  — •  -<  r- 

M)V)  S3  M  a  O)  -f 

0  0  r- 

3  3 

m  r 

f 

U  M  > 

-4 

3 

D 

aJ 

Ai 

4  4 

n  n  aj  n 

O  3  —4  4 

4  -* 

-*  -4  -4  •  •'  •-4  "C  O 

O  AJ  O  O 

0  >0  4  4 

4  4 

0  Al 

O  -O 

AJ 

p  p  n 

n  -» 

AJ 

•n 

O 

3 

M 

AJ 

O 

1 

P 

1 

4  4 

1  1 

ft  n  -t  -f 

1  t  1  1 

-4  — 1  AJ  AJ  iM 

1  1 

CM  4 

4  4  4  4  0  IO  'O 

ft  ft  ft  i'O 

AJ  IM  AJ  AJ 

AJ  AJ 

AJ  — H 

—4  —4 

- 

AJ  iM  P 

p  -1 

—1 

CvJ 

ft 

4 

4 

r-  (- 

- 

*— 4 

ft 

ft 

ft  »n 

n  ft  >n  -n 

n  0  ft  ft  n 

•n  >m 

aj  \j  aj  aj  n  ft  n 

ft  ft  ro  ■n 

ft  m  ft  ft 

n  in 

n  n 

n 

ft  ft  n 

n  aj 

M 

M 

M 

M 

M 

IS)  .M 

‘M 

AJ 

T 

C 

Y  y 

000O 

O  0  ^  C  C* 

«-  ' j 

C  0  0  0  0  0  3 

v3  O  p  : 

<  ■  C  ' '  0 

Ct  v' 

0  0 

0 

opc 

c  9 

O 

C 

c  • 

»  0  -o 

O 

4 

4 

GO  O) 

S  r-  0  -• 

i«  \  J*  3 

4  <Nl  c  j  4  CM  At  cn  «r* 

1/  <3  n  n 

—4  -*  -t  4 

t  n 

i  .4 

•  M  M 

.  M 

r-  S  0 

4  p 

0 

CO 

ft 

l> 

O  O  Os 

CO 

P 

— 4 

—4 

LA  Ut 

M  M  O 

0  0-4-40 

■n  0 

OOO  O  Vj  CM  r-4 

-4  O  O  SO  0  ;t  It 

0  n 

n  -4 

.-4  -4 

-  4 

it  it  ai 

M  X) 

D 

O 

> 

> 

—4  —4 

O' 

A- 

>- 

4 

4 

o  o 

-t  "•/)  d 

0  0  >  1>  0 

0  0 

0  0  O  O  'O  O  0 

O  O  O  A) 

ut  n  c  n 

A-  f'* 

A- 

a-  r*- 

r- 

-4  -4  Al 

AJ  ft 

•n  ft  >it 

P 

P 

it  n 

r-- 

A- 

r-> 

r- 

-4 

AJ  N  AJ  NO 

IIII 

vu  >o  r-  r-  >> 

1  1  1  1  1 

>4  -4 

-4.-4  .-4  Ul  in  4 

1  1  1  1  1  1 

; ; ; ; 

fO  in  ai  rj 

^  •-» 

--4  CJ 

AJ  AJ 

CM 

-4  -4  AJ 

1  1 

>'\1  -< 

•H 

-4 

-M 

■  -4 

CM  AJ 

t 

n 

\1 

*m 

Ifvl  frj 

M  M  U  V.l 

4^  4444  ft 

m  CO  ft  ft  ft  rt  ft 

CM  Cm  CSJ  CM  CM  AJ 

n  w  ft  *n  m  fft  ft  rti 

tt  in 

ft 

-4 

-4  -4 

*-» 

-1 

o 

o 

O  .) 

0  3  0  0 

3  3  0  3  ,3 

3  3 

3  3  0  0  0  3  3 

3  3  3  300  0  OOOOOO  3 

Q  O  3  O 

0  0 

O 

3 

O 

3  > 

3 

0 

rt 

m 

CM  »AJ 

C>J  Al  C\J  Ai 

4  4-4-4 

.L  J 

r-  sh  -C  h-  A-  A-  4 

4  0:  J  4 

0  M>  6  •  3 

4  4 

v3  »jS 

it 

n~  tv-  — * 

-<  r> 

A- 

it 

it 

Jt 

in 

m  ut 

in 

it 

o 

•n 

n  T-t  r*t  «r^  «r, 

4  <r  -4  4p 

n-  r-  4  r  cj 

1  iC'OvS  0  M3  1>  P'  n  Lt  LI  f-  i'' 

cm  oj  <4 

oD  in  cn  lit  Ut  rjD  MO  *=\J 

AJ  P 

0 

0 

i> 

r  0 

0 

O' 

X 

O 

o 

4  3  0  4 

AJ  -4 

-4  ,-Q  10  -4  O  0  -f 

r  r-  0  0 

’it  10  ft  0 

O  3 

in  4 

•>  ■  A 

4 

n  n  >n 

it  cp 

0 

3 

3  3 

0 

*3 

-* 

.  i  -  •* 

•n  0  in  wi  co  CO  co 

;\J  1T1  rji  -  j  A  p-\  ZJ 

r  — <  — «  -  4 

.4.1.4  -4 

-4  -4 

r-  •  4 

X)  n 

1^-  A-  ft 

ft  —4 

-4 

-H 

'■4 

-* 

-■* 

“I 


^4- 

w 


a 

j— < 


9 


B.  MATRIX  FUSELAGE  -  WING  REPRESENTATION 


Elastic  Matrix  -  The  coupled  matrix  provides  a  practical  method,  ver¬ 
ified  by  deflection  test  results,  for  coupling  the  vertical-longitud¬ 
inal  and  lateral-torsion  elastic  beam  properties.  The  elastic  'matrix 
includes  a  reference  axis,  an  elastic  axis  and  a  neutral  axis  which  do 
not  coincide,  thereby  coupling  the  two  directions.  Pitching,  and  yaw¬ 
ing  moments,  and  axial  load  are  transmitted  along  the  neutral  axis, 
located  at  bn  laterally  and  Cn  vertically  from  the  reference  axis. 
Vertical  and  lateral  loads,  and  rolling  moment  (fuselage  torsion)  are 
transmitted  along  the  elastic  axis  located  at  be  laterally  and  Ce 
vertically  from  the  reference  axis.  Loads  and  deflections  are  always 
referred  to  the  arbitrary  reference  axis,  at  each  end  of  the  elastic 
bay.  The  reference  axis  is  a  convenient  location  chosen  to  meet  the 
problem  requirements  for  a  particular  aircraft;  for  example,  waterline 
zero  on  the  centerplane  might  be  a  typical  reference  axis  location. 

The  coupled  elastic  matrix  is  presented  on  page  59. 

Frame  Racking  Matrix  -  In  the  deflection  test  reported  in  Reference  4, 
the  deformation  of  the  H-21  forward  fuselage  under  torsional  loading 
involved  appreciable  racking  (distortion)  of  the  fuselage  frames  and 
thus,  is  not  definable  by  the  simple  Saint  Venant  torsional  theory. 

In  addition,  it  was  determined  by  the  test  that  similar  frame  racking 
resulted  from  symmetrical  vertical  loading.  To  include  these  effects 
in  the  elastic  representation  of  the  forward  fuselage,  the  frame  rack¬ 
ing  matrix  was  written,  which  in  addition  to  normal  beam  bending  in¬ 
cludes  frame  distortion  due  to  torque  and  laterally  offset  vertical 
load. 

In  the  forward  fuselage  sections  where  frame  distortion  is  significant, 
the  conventional  torsional  representation  is  replaced  by  the  following 
system. 
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The  racking  coefficients  are  defined  from  the  respective  absolute 
deflections  at  station  n  and  n+1. 


Thus,  if  the  equations  for  vertical  and  lateral  translation  and  ro¬ 
tation  about  the  longitudinal  axis  are  added  to  the  existing  elastic 
matrix,  and  the  torsional  element  L/GJ  is  removed,  the  result  is  the 
frame  racking  matrix  shown  in  Figure  5.  It  was  found  that  an  adequate 
representation  of  H-21  forward  fuselage  coupling  required  additional 
frame  racking  from  vertical  loads.  This  was  accomplished  by  increas¬ 
ing  the  applied  moment  along  each  elastic  element  with  an  additional  moment 
equal  to  the  vertical  shear  force,  Fz  multiplied  by  a  lateral  arm,  a. 

In  the  analytical  model,  this  would  involve  a  matrix  prior  to  each  frame 
racking  matrix  increasing  the  roll  moment,  Mac  to  M©*  +  Fzd  and  then, 
after  each  elastic  element  decreasing  the  roll  moment  to  the  original 
These  matrix  operations  are  performed  in  general  terms  and  the  resulting 
matrix  illustrated  in  Figure  5.  In  the  frame  racking  matrix  as  shown, 
the  slope  is  defined  by  the  vertical  deflections  which  provide  the  best 
test  correlation. 


Coupled  Suspension  Matrix  -  The  uncoupled  fuselage  analyses,  vertical - 
longitudinal,  and  lateral-torsion  reported  in  Phase  Ila  Reference  10, 
illustrated  the  importance  of  the  engine  flexibility.  In  the  previous 
analyses  only  the  uncoupled  frequencies  were  included,  but  the  engine 
shake  test  showed  the  engine  modes  to  be  highly  coupled.  It  was  there¬ 
fore  considered  essential  that  the  suspended  engine  be  defined  by  coupled 
frequencies  and  modes.  Using  the  notation  shown  on  page  73  ,  the  mass 
motion  with  respect  to  the  fuselage  can  be  expressed  as 

e-,  =  X  e, 

Adding  the  fuselage  and  engine  motions,  the  kinetic  energy  is 

.2 


T  =?£  mi  CVe-0 


Or 


T  = 


f-  C  +  2  T  ei'1  H  y  +  X  eT  Hy) 

1*1  v  1  1  V'l  '  T-  =  |  ' 


The  total  potential  energy  of  the  system  can  be  expressed  in  terms 
of  the  effective  spring,  Kr  in ^ach  mode. 

^  ''I.  Ky 

Applying  Lagrange's  Equation  and  assuming  the  orthogonality  of  normal 
modes , 

QO  V  ^  tV}  0  0 


H 


V-l 


tn-} 


T 


Where 


M 


T 


y--\ 


mi  e> 
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Solving  the  Hr  equation  and  substituting  in  the  equation  yields, 


where 


A  1-,-j  Sj  U)z  -  Q; 

A  v„ . . 

A  1,r  z 


The  general  force  equation  derived  in  detail  in  Appendix  A,  is  expanded 
to  the  matrix  representation  for  the  coupled  engine  suspension  shown  on 
page  79. 


Uncoupled  Sprung  Mass  Matrix  -  To  provide  input  generality,  the  analysis 
includes  an  uncoupled  sprung  mass  matrix  which  is  derived  in  Appendix  A. 
Inasmuch  as  the  assumed  modes  are  uncoupled,  the  simplified  system  con¬ 
sists  of  a  suspended  mass  mounted  on  six  simple  springs  which  act  in  the 
same  directions  as  the  general  displacements,  but  are  displaced  from  the 
beam  station  by  a,  b,  and  c,  longitudinal,  lateral  and  vertical  displace¬ 
ments,  respectively. 


Using  the  amplification  factor^/  and  the  uncoupled  frequency  in  each  of 
the  fuselage  directions,  the  force  expressions  are: 


Fx  n+i  1  Fxn  +•  rru/xu^  (x  4  -  btf) 

Fy  m\  -  Fytt  +  ^2(^j"CoC  +  CL  ) 

F2  nt\  *'  Fzn  4-  no  u]2  (  z  -r  b  ^  -  &P) 

+- 

+  Q£)b  4-  Ioc^oc  U)*<* 

M^+l  -  +  vru^ui2  (x+-c£-b*)c  + 

rn Mi  (z  +  b<*  -  afi)  a  +  i?x[p  b2  $ 

Msn+l  =  (X  +  Cf  -b*)b^- 

nu^ui1  (4-  coc  +  a*)  a  * 

These  force  expressions  are  shown  in  matrix  form  as  the  uncoupled  sprung 
mass  matrix  on  page  72. 
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Mass  Matrix  -  Each  concentrated  mass  rigidly  attached  to  the  fuselage 
structure  is  included  in  the:  analytical  model  by  a  mass  matrix.  The 
mass,  with  its  associated  inertial  properties  I«x_  ,  I^  ,  and  Ia 
located  at  a,  b,  and  c,  longitudinal,  lateral  and  vertical  distances 
is  accelerated  by  karmionic  oscillation  and  produces  inertia  loads.  As 
derived  in  Appendix  A,  the  following  expressions  show  the  inertia  load 
balance  at  the  mass  station. 

I'xd+i  *'  '’xd  +  + 

^  n+i  -  *  'fifl  VO2,  ^  -j  n  -  d  +  b 

^(i+i  =  +  fn  ui  (  £x\ +  b  <x  n  -  Q_  ft) 

M<xnH  =  -  rntu Sz  (y»»-  cocn  +  a^c  + 

nruii2- (^-n -h  boc^  -  a^n)b  ^  I  c*n 
Mg  o +\  =  f\/| $ n  +■  nn u32  b  C  ~ 

who*  (z„  botn  -aWdfi, 

ntl  =  4’ 

+  TTi'b  (kLjn"C,cXn+  +•  X  ^  VAj2  ^  yi 

These  fuselage  load  equations^  together  with  the  unchanged  fuselage 
deflections  at  the  station,  are  expressed  in  matrix  form  to  produce 
the  mass  matrix  shown  on  page  61-  l 

Ground  Spring  Matrix  -  External  spring  restraints  such  as  used  in  the 
ground  shake  test  for  suspension  of  the  aircraft  from  the  roof  trusses 
of  a  hangar  building  are  simulated  in  an  analytical  model  by  the  inclu¬ 
sion  of  a  ground  spring  matrix  at  each  support  point.  If  a  beam  station 
is  restrained  by  springs,  additional  loads  are  induced  on  the  structure 
in  proportion  to  the  motion  in  the  direction  of  each  spring.  In  matrix 
form  the  force  expression  shown  in  Appendix  A,  combined  with  the  constant 
fuselage  directions,  results  in  the  ground  spring  matrix  illustrated  on 
page  63. 

Concentrated  Spring  Matrix  -  Concentrated  springs  between  adjacent  beam 
sections  are  necessary  in  the  analytical  representation  for  simulating 
localized  flexibility  of  the  structure,  such  as  the  attachment  between 
the  rotor  transmission  and  the  fuselage  structure.  Since  no  inertia 
loads  are  involved,  the  only  changes  across  the  joint  are  deflections 
resulting  from  the  action  of  the  loads  on  the  spring.  Therefore,  the 
loads  can  be  equated  to  the  spring  rate  multiplied  by  the  deflection 
change  as  shown  below.  r_ 

~  Ji'  ^  X  YHl  ^  =  Xyi  + 

yn  -  +  (  ^r\-H  “Mr\\  j  n+i  -  ^ 


Emil 

kT 
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f"Hn  =  +■  (  Z’n-vi  -Zn]  a  ^n  +  i  *  +  ^~2l[y/Kz. 

M  *n  =  F  Koc  +  ^  OC  r> 4- ,  =  n  +■  Mocn^ 

Mj^T)  s  f  ^j?  (  £vhi  "  $n)  i  i^n+i  =  4  Mjw 

Mxn  r  +  K*  (  Xri+-t  ‘  Xn)  )  T 

The  concentrated  spring  matrix  on  page  65  is  obtained  using  these  dis¬ 
placement  expressions  together  with  the  condition  of  no  load  change 
across  the  joint. 


Vertical  Bend  Matrix  -  A  general  matrix  program  for  elastic  structures 
must  permit  a  change  in  direction  of  the  elastic  axis  such  as  used  for 
a  "U"  beam  representation  of  the  H-21  fuselage.  In  progressing  from 
station  n  to  n+1  the  loads  and  deflections  are  oriented  from  the  orig¬ 
inal  position  to  a  direction  corresponding  to  the  new  elastic  axis 
location.  The  general  procedure  is  illustrated  below  for  the  longitud- 


inal  Fx  and 

vertical  Fz  forces. 

E.  ft. 

,  Fzn 

- *-  Fxri 

Fxn  +  t  = 

Fxvi 

Cos  $  -  F?n 

Bin  $ 

Fzn+t  = 

Fxn 

5  i  n  B  +  Fzn 

Cos 

r>+l 


Using  the  same  procedure  for  the  additional  loads  and  deflections  and 
considering  only  rotation  in  the  vertical  plane,  the  vertical  bend 
matrix  on  page  69  is  obtained. 


Shift  Matrix  -  Inclusion  of  the  engine  suspension  in  the  matrix  program 
requires  a  transformation  of  loads  and  deflections  from  station  521,  the 
point  of  fuselage  attachment,  to  the  engine  C.G.  at  which  the  engine 
natural  modes  are  defined.  Considering  a  station  (rrt-1)  located  a,  longi¬ 
tudinally,  b,  laterally  and  c,  vertically  from  the  original  station,  the 
transformation  equations  shown  below  define  the  loads  and  deflections  at 
the  new  position. 

Fxn-H  *  Fxn  -  n+.i  =  ^ 

F y  rvt-i  =  Fy  n  $  y\  v  t  *  $  h 

Fan+i  *  ^vut  = 

K/)o<  p  +. |  =  M -  bo  Fg  ^  +  Co  Ftjn  y  ^ nw  =  ^ n  -  Fo  Co 

n4M  =  Mpn  ~  Co  Fxn  ■*"  ^  ^-jn+i  -  ^ 

K/|  =■  M  Xn  +  bo  Fxn  ”  Co  bo  °^r\  ”  Q-* 

In  matrix  form,  these  expressions  are  shown  as  the  shift  matrix  on 
page  67, 


Forced  Matrix  -  In  addition  to  the  beam  matrices  necessary  for  a 
natural  mode  analysis,  the  forced  solution  requires  one  additional 
matrix  which  permits  a  forcing  function  to  be  inserted  at  any 
position  on  the  structure.  Considering  that  external  loads  are 
applied  to  a  beam  station  n,  the  loads  acting  at  station  n+1  are 
shown  below  as  derived  in  Appendix  A. 


c 

!  x  n-H  " 

Fxn 

4.  cr 

•  i  X 

Fm  m-H  = 

Fcjn 

+  Fy 

Fztui  = 

FZp 

+  F* 

M<*  n 

-  C  Rj  ■ 

*  bFz 

+ 

Mpn+i  5 

•+  C  Fx 

-  a 

Fb 

M  l  n+i  = 

-  bFx 

+  a 

F„ 

k. 

vu_ 

+  M  p 

*  Mx 


In  the  forced  frequency  solution,  the  forcing  functions  are  con¬ 
stant  loads  which  are  added  to  the  inertia  and  elastic  loads.  To 
add  the  constant  values,  the  analysis  is  expanded  to  a  13  x  13  matrix 
and  thus,  permits  the  addition  of  a  constant  at  any  beam  station.  Ex 
pressing  the  forcing  function  as 

F  s  Fo  +  F,  co5i|jn  +  Fz  5\MVUr» 

Where,  r\  -  Harmonic  Azimuth 

The  forced  matrix  can  be  written  as  shown  on  page  79. 
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f - 

Right  Wing 

l  tVi  *  --  — •  ''  1 

Effective  Flap  Hinge- 


£  Fuselage 


Considering  the  fuselage-wing  schematic  shown  above,  the  right  wing  ele¬ 
mental  mass  motion  as  derived  in  Appendix  A  is, 


=  x  -  b/tf 
=  y  +  a^  & 

=  z  =  bjf*-  atP>+  zt  -  a1^  ± 
where,  Z^,  o<  ^  Wing  Displacements 

x,  y,  z,  c<  , ^  Fuselage  Displacements 
Summing  along  the  wing,  the  total  kinetic  of  the  right  wing  is 

"(x-b^f 
+  a,  *) 

+  (z  +  "  a-J  *  -  Qi  A)2 

zK  c*  *  -,t  -  ift  ( w  *  i«  (v  *  < 


Applying  Lagrange's  Equations,  the  right  wing  equations  of  motion  are: 

x  C?miVx  -  Fx 

(I  mOy  +  (  ?  mi  at)  Y  -  fs 

(X  m j)z  +  (  X  rri;  bi)£  +  (-^ rrii  a-,)  } 

Fz 


y 

E 


_  X  ( ^  *  let)  «•  •  (2-  m,  bi)H  +(-^.  rri,  Qv  \\)$ 
( m-,  al  +  i  )J  0  +■  (-y  rri-,  a-t)  z  +  (-Z  m-t  q;  b,) «. 


Hs  +  Ks  Hs 


V  ^  /  (S)  aCS)-]  Oo 

Lym>6  -aj Je  + 

V  -r 

o  C 

=0 

0  o 


Assuming  harmonic  motion  and  solving  for  Hg 

i_j  =  £ ^  ■+-  c\  i-  q^s  g  )  {Jj2 


1  »s 


^I3-^IS 
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where, 


O^s 


Oi. 


I. 


O'. 


7S 


‘  ms  a-t  (2;s>  -  +  X  I?i  & 

■  ^ms^czr-  q.o 

1  ?  (?“-  QiCf  +  ?  ^>2+?  I=,i  oi 

•  £  m-,  0*iM- *$“) 


s)2 


Substituting  Hg  in  the  equations  of  fuselage  motion  the  right  wing 
inertial  loads  can  be  expressed  in  matrix  form  as  shown  on  page  '  2 . 

In  a  similar  manner,  the  left  wing  inertial  loads  are  obtained  and 
added  to  the  right  wing  resulting  in  the  expressions  shown  on  page  93. 
Expanding  to  the  general  matrix  form,  the  effective  wing  matrix  is 
illustrated  on  page  94. 
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SECTION  III 


WING-FUSELAGE  FORCED  RESPONSE 


A.  GENERAL 

Forced  response  of  the  H-2I  belicoPter  equipped  with  floating  fuel  wings 
is  calculated  using  measured  hub  load  data  from  Reference  7.  These  loads 
were  measured  on  an  H-21  helicopter  and  are  applied  to  the  analytical  rep¬ 
resentation  of  the  H-21  fuselage-wing  combination  so  as  to  provide  reason¬ 
able  prediction  of  the  forced  response.  Inclusion  of  the  floating  fuel 
wings,  unlike  the  configuration  used  in  the  load  measurement  program  intro¬ 
duces  some  error  in  regard  to  applied  loads.  However,  this  effect  becomes 
significant  only  if  a  resonant  condition  exists  with  the  forcing  frequency, 
which  could  be  evident  in  the  undamped  response. 

As  in  the  three-bladed  H-21  helicopter  without  fuel  wrings,  third  harmonic 
vibration  levels  are  usually  a  dominant  portion  of  the  vibration  environ¬ 
ment;  therefore,  calculations  are  performed  using  only  the  fixed  system 
third  harmonic  loads.  Hub  load  data  in  Reference  7  was  determined  in  two 
forms;  fixed  system  shaft  loads  which  include  rotor  inertia  and  damping 
effects,  and  fixed  system  aerodynamic  loads  which  are  the  shaft  loads  ad¬ 
justed  for  inertia  loads  resulting  from  fuselage  induced  rotor  hub-motion. 
In  this  analysis,  the  rotor  shaft  loads  are  used,  thereby  eliminating  the 
complexity  of  the  rotor  system  dynamics  in  the  analytical  representation. 

Forced  response  of  the  H-21  fuselage-wing  combination  is  investigated  at 
maximum  cruise  speed  of  90  knots  for  three  rotor  speeds  within  the  opera¬ 
ting  rotor  speed  band.  The  analysis  considers  third  harmonic  loads 
corresponding  to  rotor  speeds  of  250,  ?58  and  268  RPM.  In  addition,  the 
two  extremes  of  the  wing  are  considered,  which  are  OX  fuel  and  100X  fuel 
corresponding  to  approximate  weights  of  1000  lbs.  per  wing  and  8000  lbs. 
per  wing  respectively.  Considering  the  extreme  weight  limits  of  the  wing, 
the  forced  modes  provide  reasonable  approximations  of  the  third  harmonic 
vibration  environment  which  the  H-21  helicopter-wing  combination  would 
experience  under  normal  flight  conditions. 

B.  ANALYTICAL  MODEL 

Fuselage  -  Previous  theoretical  studies  performed  under  Reference  8,  9, 
and  10  showed  the  requirement  for  and  subsequent  development  of  an 
analytical  representation  of  the  H-21  fuselage  which  included  vertical- 
lateral  coupling.  Figure  6  illustrates  the  final  representation  of  the 
H-21  fuselage.  In  the  forward  fuselage,  each  elastic  section  is  typified 
by  a  frame  racking  matrix  which  in  addition  to  normal  beam  bending  in¬ 
cludes  frame  distortion.  Figure  7  presents  the  complex  elastic  properties 
for  each  forward  fuselage  section.  The  adequacy  of  the  forward  fuselage 
representation  is  demonstrated  by  Figure  8  which  compared  the  calculated 
forward  fuselage  deflection  under  a  3600  lb.  vertical  load  to  the  deflec¬ 
tion  data  obtained  in  the  fuselage  load  test. 
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SECTION  III 


WING -FUSELAGE  FORCED  RESPONSE 


A,  GENERAL 

Forced  response  of  the  H-21  helicopter  equipped  with  floating  fuel  wings 
is  calculated  using  measured  hub  load  data  from  Reference  7.  These  loads 
were  measured  on  an  H-21  helicopter  and  are  applied  to  the  analytical  rep¬ 
resentation  of  the  H-21  fuselage-wing  combination  so  as  to  provide  reason¬ 
able  prediction  of  the  forced  response.  Inclusion  of  the  floating  fuel 
wings,  unlike  the  configuration  used  in  the  load  measurement  program  intro 
duces  some  error  in  regard  to  applied  loads.  However,  this  effect  becomes 
significant  only  if  a  resonant  condition  exists  with  the  forcing  frequency 
which  could  be  evident  in  the  undamped  response. 

As  in  the  three-bladed  H-21  helicopter  without  fuel  wjngs,  third  harmonic 
vibration  levels  are  usually  a  dominant  portion  of  the  vibration  environ¬ 
ment;  therefore,  calculations  are  performed  using  only  the  fixed  system 
third  harmonic  loads.  Hub  load  data  in  Reference  7  was  determined  in  two 
forms;  fixed  system  shaft  loads  which  include  rotor  inertia  and  damping 
effects,  and  fixed  system  aerodynamic  loads  which  are  the  shaft  loads  ad¬ 
justed  for  inertia  loads  resulting  from  fuselage  induced  rotor  hub-motion. 
In  this  analysis,  the  rotor  shaft  loads  are  used,  thereby  eliminating  the 
complexity  of  the  rotor  system  dynamics  in  the  analytical  representation. 

Forced  response  of  the  H-21  fuselage-wing  combination  is  investigated  at 
maximum  cruise  speed  of  90  knots  for  three  rotor  speeds  within  the  opera¬ 
ting  rotor  speed  band.  The  analysis  considers  third  harmonic  loads 
corresponding  to  rotor  speeds  of  250,  ?58  and  268  RPM.  In  addition,  the 
two  extremes  of  the  wing  are  considered,  which  are  0%  fuel  and  100%  fuel 
corresponding  to  approximate  weights  of  1000  lbs.  per  wing  and  8000  lbs. 
per  wing  respectively.  Considering  the  extreme  weight  limits  of  the  wing, 
the  forced  modes  provide  reasonable  approximations  of  the  third  harmonic 
vibration  environment  which  the  H-21  helicopter-wing  combination  would 
experience  under  normal  flight  conditions. 

B.  ANALYTICAL  MODEL 


Fuselage  -  Previous  theoretical  studies  performed  under  Reference  8,  9, 
and  10  showed  the  requirement  for  and  subsequent  development  of  an 
analytical  representation  of  the  H-21  fuselage  which  included  vertical- 
lateral  coupling.  Figure  6  illustrates  the  final  representation  of  the 
H-21  fuselage.  In  the  forward  fuselage,  each  elastic  section  is  typified 
by  a  frame  racking  matrix  which  in  addition  to  normal  beam  bending  in¬ 
cludes  frame  distortion.  Figure  7  presents  the  complex  elastic  properties 
for  each  forward  fuselage  section.  The  adequacy  of  the  forward  fuselage 
representation  is  demonstrated  by  Figure  8  which  compared  the  calculated 
forward  fuselage  deflection  under  a  3600  lb.  vertical  load  to  the  deflec¬ 
tion  data  obtained  in  the  fuselage  load  test. 
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FIGURE  7{cONT)H-2i  FORWARD  FUSELAGE  ELASTIC  PROPERTIES 


Figure  9  shows  the  elastic  sections  used  for  representing  the  aft  fuselage 
Weight  distribution  for  the  dynamic  model  of  the  fuselage  is  shown  in  Fig¬ 
ure  10,  together  with  the  fuselage  reference  axis. 

As  a  relatively  important  dynamic  element,  the  fuselage  model  includes  the 
engine  as  an  elastically  supported  mass.  An  engine  shake  test,  reported 
in  Reference  8,  was  used  to  obtain  the  coupled  modes  and  frequencies  illus¬ 
trated  as  Figure  11,  for  use  in  the  analytical  fuselage  representation 

Floating  Fuel  Wings  -  Each  floating  fuel  wing  is  attached  to  the  helicopter 
through  a  skewed  hinge  so  as  to  eliminate  the  bending  moment  applied  at  the 
fuselage  by  a  conventional  wing.  The  analytical  representation  of  each  wing 
includes  rigid  body  motion  to  simulate  rotation  about  the  hinge  and  all 
natural  modes  up  to  and  including  the  natural  in  the  vicinity  of  3  iTt excita¬ 
tion.  Using  the  stiffness  and  mass  properties  shown  in  Figure  12,  calcula¬ 
tions  are  performed  to  obtain  the  wing  natural  frequencies  and  modes.  Figure 

13  illustrates  the  wing  modes  for  the  07.  fuel  configuration  similarly.  Figure 

14  presents  the  wing  modes  for  the  1007.  fuel  configuration.  The  illustrated 
wing  modes  show  only  the  semi-span  shapes,  but  as  a  result  of  the  coupled 
vertical-lateral  motions,  both  symmetrical  and  anti-symmetrical  modes  exist 
at  the  indicated  natural  frequencies.  Representing  each  of  the  uncoupled 
wing  modes  as  generalized  coordinates,  the  dynamic  wing  loads  are  transmitted 
to  the  fuselage  in  the  coupled  wing-fuselage  analysis  using  the  effective 
wing  matrix. 

C.  FORCED  RESPONSE 


Forced  Mode  Solution  -  Using  the  analytical  representation  of  the  H-21  heli¬ 
copter  equipped  with  floating  fuel  wings,  the  forced  responses  are  obtained 
from  the  coupled  matrix  analysis  programmed  for  the  Univac  1103A  digital 
computer.  A  force  matrix  at  each  rotor  applies  the  measured  loads,  and  the 
matrix  array  is  collapsed  using  the  forcing  frequency  to  apply  acceleration 
to  each  mass  and  inertia  item,  resulting  in  an  applied  load  and  moment  dis¬ 
tribution.  Applying  the  free-free  boundary,  the  fuselage  matrix  collapses 
to  a  6  x  6  boundary  determinant  which  is  solved  for  the  forward  rotor  de¬ 
flection.  For  each  forcing  frequency  and  the  corresponding  measured  loads, 
the  forced  response  is  obtained  by  a  detailed  listing  of  intermediate  results 
of  each  matrix  multiplication  starting  with  the  forward  rotor  boundaries  and 
continuing  along  the  reference  axis  to  the  aft  rotor. 

Forced  Mode  Calculations  -  Using  the  analytical  representation  of  the  H-21 
equipped  with  floating  fuel  wings  and  the  measured  third  harmonic  loads 
illustrated  in  Figure  15,  matrix  calculations  are  performed,  and  the  results 
presented  in  Figures  16  to  18,  and  20  to  22,  as  pIhiwh  below. 
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FIGURE  10  -  H-21  REFERENCE  AXIS  LOCATION  AND  WEIGHT  DISTRIBUTION 


Modes  where  angles  dominate  have  those  angles  normalized  to  0.04  radians; 
0.04  radians  “  1"  deflection  at  25"  radius.  Where  displacements  dominate 
the  displacements  are  normalized  to  one  inch. 


Figure  11.  Normalized  Engine  Modes,  H-21C-96  Engine  Shake  Test 
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FIGURE  13  WING,  FLAP  BENDING  AND  TORSION  MODES  - 


Third  Harmonic  Forced  Response  -  0‘/.  Wing  Fuel 


For  the  07.  wing  fuel  case,  time  histories  of  third  harmonic  forced  re¬ 
sponse  are  illustrated  in  Figures  16,  17  and  18  during  an  RPM  sweep  at 
90  knots.  Each  figure  presents  the  instantaneous  fuselage  and  engine 
amplitudes  at  four  discrete  time  positions  during  a  half  cycle  of  third 
harmonic  oscillation.  The  harmonic  azimuth  positions  shown  are  at  0,  45, 

90  and  135  degrees.  The  corresponding  azimuth  positions  between  180  and 
360  degrees  repeat  the  original  response  in  the  opposite  direction. 

As  shown  in  the  vertical  bending  plots,  the  forward  and  aft  pylon  longi¬ 
tudinal  deflections  are  normally  in-phase,  i.e.,  motions  occur  simultane¬ 
ously  in  the  same  direction.  The  longitudinal  magnitudes  increase  with 
RPM,  reaching  peaks  of  0.08  inches  at  Che  forward  totor  and  0.12  inches 
at  the  aft  rotor  at  268  RPM.  In.  the  vertical  direction,  the  peak  response 
occurs  for  fuselage  station  603  at  a  rotor  speed  of  258  RPM  with  a  180 
degree  phase  shift  from  the  amplitude  at  250  RPM  and  then,  decreases  to 
a  negligible  vertical  response  at  268  RPM. 

Laterally,  at  250  RPM,  the  fuselage  response  is  small.  However,  when  the 
rotor  speed  is  increased  to  258  RPM,  the  fuselage  responds  predominantly 
in  rigid  yaw  which  progresses  to  the  first  lateral  bending  mode  at  268  RPM. 
Maximum  lateral  amplitudes  of  0.09  and  0.10  inches  occur  near  Station  100 
and  at  the  maximum  rotor  speed.  Generally,  the  torsional  response  of  the 
fuselage  is  small,  although  the  aft  fuselage  shows  some  frequency  sensitive 
response  that  increases  with  rotor  speed. 

Also,  during  the  RPM  sweep  at  90  knots,  the  engine  exhibits  significant 
changes  in  amplitude.  In  the  lateral  direction,  the  engine  response  is 
negligible  at  250  RPM,  but  increases  with  RPM,  generally  out-of-phase  with 
the  fuselage  lateral  motion,  to  a  maximum  of  0.08  inch  at  268  RPM.  Yaw 
motion  of  the  engine  shows  an  apparent  phase  reversal  between  250  and  258 
RPM,  but  relative  to  the  aft  fuselage  vertical  motion  the  phase  is  constant. 
Similarly,  the  vertical  motion  is  out-of-phase  with  the  fuselage  and  follows 
the  vertical  amplitude  trend  of  the  aft  fuselage  with  a  peak  response  of 
0.20  inch  at  258  RPM.  It  is  of  interest  to  note  that  at  258  RPM  the  longi¬ 
tudinal  engine  motion  is  in-phase  with  the  pylons,  but  out-of-phase  with 
the  pylons  at  the  higher  and  lower  operating  speeds. 

The  fuselage  responses  shown  are  influenced  by  the  natural  frequencies 
of  the  wing-fuselage  combination,  illustrated  for  the  fuselage  without 
wings  in  the  matrix  residual  curve  of  Figure  19.  Inclusion  of  the  07. 
fuel  wing  frequency  at  10.39  CPS  would  influence  the  residual  curve  by 
producing  an  additional  crossing  corresponding  to  the  wing  mode  and  shift¬ 
ing  the  fuselage  frequencies  to  include  wing  coupling.  However,  as  shown 
in  the  mode  plots  the  wing  does  not  significantly  restrain  the  motion  of 
the  forward  fuselage  because  of  its  light  weight  relative  to  the  fuselage 
alone  Therefore,  it  is  reasonable  to  consider  that  the  wing  mode  does 
not  appreciably  influence  the  fuselage  forced  response  in  the  0%  fuel 
configuration. 


34 


THIRD  HARMONIC  FORCED  RESPONSE 


0%  WING  FUEL,  90  KNOTS,  268  RPM 
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FIGURE  18 


FIGURE  19  MATRIX  RESIDUAL  PLOT 


Third  Harmonic  Forced  Response  -  100%  Wing  Fuel 


Figures  20,  21  and  22  illustrate  the  third  harmonic  forced  response 
during  an  RPM  sweep  at  a  cruise  speed  of  90  knots.  Presented  in  the 
same  format  as  the  previous  /response  plots,  the  instantaneous  fuselage 
and  engine  amplitudes  are  shown  for  harmonic  azimuth  positions  of  0, 

45,  90,  and  135  degrees. 

At  250  RPM  in  Figure  20,  longitudinal  pylons  in  phase  response  is  seen 
to  occur,  but  the  amplitude  relative  to  lateral  fuselage  motion  is  much 
smaller  than  it  was  during  the  0%  fuel  copdition.  Lateral  motion  is 
very  large,  with  lateral  forward  fuselage  motions  of  +  0.5  inches  at 
3.0. ,  and  correspondingly  large  torsional  amplitudes.  Note  that  the 
scales  in  Figure  20  differ  from  the  other  figures  because  of  the  large 
lateral  amplitudes.  The  fully  fueled  wing  is  not  easily  excited  because 
of  its  large  mass  so  that  the  wing  point  at  Station  259  is  a  node,  with 
the  fuselage  bending  about  this  point  as  though  it  were  a  pivot  magnify¬ 
ing  the  cockpit  floor  lateral  motion.  This  large  response  is  indicative 
of  a  natural  frequency  of  the  coupled  fuselage-wing  system. 

Natural  modes  of  the  helicopter  without  wings  are  represented  by  the  zero 
crossings  in  the  Figure  19  residual  curve.  Introduction  of  the  wing 
natural  modes  of  Figure  14  as  asymptotes  would  modify  the  residual  plot 
and  produce  additional  coupled  natural  frequencies.  Inasmuch  as  the 
residual  curve  was  calculated  for  a  gross  weight  of  11,100  lb  whereas 
the  helicopter  wing  system  with  100%  fuel  has  a  gross  weight  of  27,100  lb 
the  natural  modes  shown  on  the  residual  curve  provide  no  definite  verifi¬ 
cation  of  the  predicted  resonance. 


However,  the  previous  matrix  calculations  without  the  fuel  wings  as  rep¬ 
resented  on  the  residual  curve  show  a  lateral  fuselage  mode  at  18.9  CPS 
with  maximum  amplitude  in  the  vicinity  of  the  wing  attachment.  Using  the 
Jrequencjr  correction  as  the  square  root  of  the  gross  weights, 

x  18.9,  the  corresponding  natural  frequency  of  the  fuselage¬ 


(11.100)1% 

(27,100) 


wing  system  is  12.1  CPS.  Therefore,  it  is  probable  that  this  natural  mode 
near  3d  results  from  the  18.9  CPS  lateral  fuselage  mode 


Increasing  RPM  reduces  the  lateral  response  amplitude  as  shown  in  the 
Figure  21  plot  at  258  normal  RPM.  Vertical  amplitudes  are  now  of  about 
the  same  order  as  the  laterals,  with  rotor  pylon  peaks  near  t  o.l  inches 
The  scales  in  Figure  21  are  enlarged  back  to  their  previous  values,  but 
the  plots  of  course  represent  smaller  amplitudes.  The  highest  RPM  con¬ 
sidered,  268  RPM  in  Figure  22  exhibits  similar  vertical  amplitudes,  but 
lateral-torsion  amplitudes  are  further  reduced. 
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THIRD  HARMONIC  FORCED  RESPONSE 
100%  WING  FUEL,  90  KNOTS,  268  RPM 


1 1  {-] :  ! ; 1  '  i  :  '•4-  *  U_L.iL-._l  _il . 

TTh-  | ;  o,  i ;  •  •  • :  1 

\.\  Liill;  f.'UHiT.*  *  ,0  Till,,  «* 

. aiMi  wtomimm 


AZIMUTH  - 


f *0  AVION 
(LOOKING  AFT) 


0.1  -0.1 


AFT  FVIOH 
(LOOKING  AFT) 


•  '  I  AOTOA  Hf  AD 
M  ''  TOA  OF  «MSN 


»TA.  TOO  200  300  400 


m 


azimuth  -  1350  i-. 

WmSflfi  liflw 

mHMKI 

wimfk mMllJi  ■  i :  ■.  \  •  r  -  *  ^  -  -  . 


■'Uitiii 


■si 


0.1  -0.1 
F«B  AVION 
(LOOKING  AFT) 


0.1  -u* 
AFT  AYLON 
(LOOKING  AFT) 


0.1  -0.1 
F*D  AVION 
(LOOKING  AFT) 


2»-1 


FIGURE  22 


Third  Harmonic  Cockpit  Floor  Motion 


Calculated  cockpit  floor  responses  for  the  standard  fuselage  without 
wings  are  compared  in  Figure  23  with  measured  flight  amplitudes.  This 
figure,  taken  from  Reference  8  ,  presents  calculated  amplitudes  as 
solid  lines  and  measured  amplitudes  as  dashed  lines.  Vertical  calcu¬ 
lated  amplitudes  compare  very  well  with  measured  amplitudes.  Lateral 
amplitude  is  satisfactory  at  240,  250  and  278  RPH,  but  compare  poorly 
at  260  and  270  RPM. 

Figure  24  presents  similar  cockpit  floor  calculated  amplitudes  for  the 
helicopter  with  wings.  In  the  vertical  direction  with  0%  wing  fuel, 
the  amplitude  at  each  operating  speed  is  nearly  identical  to  those  cal¬ 
culated  for  the  helicopter  without  wings.  Laterally  for  the  same  con¬ 
figuration,  at  250  RPM,  the  responses  are  again  nearly  identical  with 
an  amplitude  of  .002  inch,  but  with  increasing  rotor  speed  the  response 
of  the  helicopter  with  wings  increases  above  that  of  the  standard  heli¬ 
copter  Calculations  show  a  peak  response  of  0.10  inch  at  268  RPM 
approximately  five  times  the  response  of  the  helicopter  without  float¬ 
ing  fuel  wings. 

For  the  1007.  fuel  configuration,  the  maximum  undamped  cockpit  floor 
response  occurs  in  the  lateral  direction  at  250  RPM.  From  a  calculated 
peak  amplitude  of  0J64  inch,  the  lateral  response  decreases  as  the  rotor 
speed  increases  reaching  0.17  inch  at  258  RPM  and  then,  at  268  RPM  a 
minimum  value  of  0.005  inch  close  to  the  calculated  response  of  the 
helicopter  without  wings.  Corresponding  to  the  near  resonant  condition 
in  the  lateral  direction,  the  .035  inch  vertical  peak  occur  at  250  RPM. 

As  the  rotor  speed  is  increased,  the  vertical  response  appears  similar 
to  that  calculated  for  the  helicopter  without  wings. 

The  very  large  amplitudes  calculated  here  at  low  RPM  indicate  the  rapid 
buildup  of  amplitude  for  an  undamped  resonance.  From  past  experience, 
although  large  vibration  levels  can  occur,  amplitude  increases  as  large 
as  that  in  Figure  24  are  unlikely  in  the  real  damped  case.  The  results 
do  indicate  that  the  heavy  fueled  wing  will  tend  to  induce  higher  than 
normal  fuselage  vibration  levels  during' the  early  part  of  a  ferry  mission, 
until  the  fuel  quantity  is  somewhat  reduced.  It  is  also  clear  that  there 
is  a  strong  tendency  toward  higher  levels  with  reduced  rotor  RPM,  but 
from  performance  studies  the  RPM  will  tend  to  be  closer  to  or  slightly 
above  normal  for  best  range. 
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FIGURE  24  COCKPIT  FLOOR  RESPONSE  -  90  KNOT,  RPM  SWEEP 
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SECTION  IV 


CONCLUSIONS 


The  H-2I  helicopter  equipped  with  floating  fuel  tanks  for  ferry  range 
extension  has  been  investigated  for  fuselage  vibratory  response  char¬ 
acteristics  . 

Third  harmonic  response  of  the  H-21  equipped  with  floating  fuel  tanks 
was  investigated  for  the  two  extremes  of  wing  loading,  0%  and  100%  wing 
fuel,  and  compared  to  the  cockpit  floor  vibration  of  the  helicopter  with¬ 
out  fuel  wings.  Calculated  using  measured  third  harmonic  loads,  the 
forced  modes  corresponding  to  07.  wing  fuel  are  strongly  coupled  lateral- 
vertical  modes  with  significant  fuselage  motion.  Cockpit  floor  motion 
compares  favorably  in  the  vertical  direction  with  that  calculated  for  the 
helicopter  without  fuel  wings,  but  laterally,  the  fuselage  response  in¬ 
creases  with  rotor  speed  to  a  maximum  of  0.10  inch  at  268  RPM.  Therefore, 
the  acceptable  cockpit  3 JTL  vibration  environment  associated  with  the  heli¬ 
copter  without  wings  is  maintained  during  the  lower  rotor  speeds  when  the 
fuel  tanks  are  nearly  empty.  At  the  higher  rotor  speed,  the  calculated 
3  Ci.  vibration  is  at  a  high  level,  but  damping  should  reduce  these  ampli¬ 
tudes  to  a  nearly  satisfactory  level. 

In  the  100%  fuel  configuration,  the  calculated  30.  vibration  level 
appears  high  at  the  lower  rotor  speeds,  but  satisfactory  at  268  RPM, 
Fuselage  mode  plots  show  the  same  trends  with  RPM  and  further,  indicate 
the  presence  of  a  lateral  fuselage  mode  in  the  vicinity  of,  but  below 
the  3 n.  excitation  at  250  RPM.  Laterally,  at  250  RPM  the  cockpit  floor 
motion  peaks  at  0.76  inch  reflecting  the  undamped  natural  mode,  and  then 
decreases  rapidly  to  an  amplitude  of  .005  inch  at  268  RPM.  Reflecting 
the  near  lateral  resonance,  the  vertical  cockpit  floor  motion  at  250  RPM 
shows  an  undamped  response  of  .035  inch.  However,  above  258  RPM,  the 
third  harmonic  vertical  vibration  level  is  nearly  equal  to,  or  less  than 
those  calculated  for  the  helicopter  without  fuel  wings.  Therefore,  with 
full  fuel  tanks,  the  undamped  calculations  show  a  satisfactory  vibration 
environment  would  exist  in  the  aircraft  at  the  maximum  rotor  speed.  At 
the  lower  rotor  speeds,  the  calculated  amplitudes  are  high  reflecting 
the  undamped  lateral  resonance.  Past  experience  has  shown  that  ampli¬ 
tude  increases  as  large  as  those  calculated  are  very  unlikely  in  the 
real  damped  case. 

It  is  recommended  that  prior  to  flight,  a  simplified  ground  shake  test 
be  performed  on  the  H-21  fuselage-wing  combination  to  substantiate  the 
theoretical  calculations.  In  addition  to  providing  a  check  on  the  cal¬ 
culated  resonance  below  the  3.0.  excitation  at  250  RPM  with  100%  fuel, 
the  test  would  provide  reasonable  prediction  of  the  fuel  and  structural 
damping.  Further,  the  proposed  shake  test  may  provide  advantageous 
partial  fuel  arrangements  and  a  sequence  for  emptying  the  fuel  tanks 
which  by  controlling  the  important  natural  frequencies  and  increasing 
fuel  damping  would  improve  the  vibration  characteristics. 
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APPENDIX  A 


COUPLED  MODE  ASSOCIATED  MATRIX  DERIVATION 


i 

COUPLED  MODE  ASSOCIATED  MATRIX  DERIVATION 
Associated  Matrix  Procedure 

For  fuselage  vertical-lateral  bending  nodes,  six  degrees  of  freedom  are  described  by 
the  matrix  terms;  longitudinal,  lateral  and  vertical  displacement,  roll,  pitch  and 
yaw  rotation.  The  structure  is  separated  into  lumped  parameter  form  and  a  matrix 
representation  prepared  for  each  property.  The  matrices  are  assembled  into  an 
array  to  simulate  a  progression  from  one  end  of  the  fuselage  beam  to  the  other. 

Starting  with  a  set  of  load  and  deflection  boundary  conditions  at  one  end  of  the 
beam,  the  top  of  the  forward  rotor;  successive  multiplication  of  numerical  matrices 
is  performed  to  reach  the  other  end  of  the  beam.  For  this  operation,  a  trial 
frequency, oo  ,  is  used,  which,  by  the  conventional  harmonic  motion  assumption, 
applies  an  acceleration  to  each  mass  and  inertia  item,  resulting  in  an  applied 
load  and  moment  distribution  on  the  beam.  At  the  other  end  of  the  beam,  a  second 
set  of  boundary  conditions  are  enforced;  by  successive  trials,  frequencies  are 
found  which  satisfy  the  boundaries  and  are,  therefore,  the  natural  frequencies 
being  sought.  Each  trial  matrix  multiplication  produces  a  non-zero  residual;  when 
this  value  is  at  or  very  near  zero,  the  natural  frequency  is  established. 

When  the  natural  frequency,^  n ,  is  determined,  one  further  matrix  multiplication 
of  the  system  is  required;  thiscon  value  and  an  assumed  unit  deflection  at  one 
boundary  is  used  to  obtain  a  detailed  listing  of  the  intermediate  results  of  each 
matrix  multiplication,  thus  providing  the  natural  mode  shapes  corresponding  to0on> 

This  listing  contains  displacements,  rotatiews,  forces  and  moments  at  each  station 
along  the  structure.  In  practice,  the  extensive  matrix  multiplications  are  performed 
on  a  digital  computer  and  lead  to  rapid  solution. 
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Elastic  Matrix 


The  general  elastic  matrix  for  a  weightless  beam  element  undergoing  deflections 
in  the  vertical  plane  was  derived  in  Appendix  A,  Reference  10 
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The  elastic  matrix  for  a  weightless  beam  element  undergoing  deflections  in  the 
lateral  plane  was  also  derived  in  Appendix  A,  Reference  10 
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The  effect  of  the  axial  force  <3  in  the  above  matrices  will  be  neglected  in  the 
formation  of  the  coupled  elastic  matrix  because  experience  has  proved  it  to  have 
little  effect.  The  first  step  in  forming  the  coupled  matrix  is  the  writing  of  a 
twelfth  order  matrix  containing  both  the  above  uncoupled  components  in  their  proper 
diagonal  positions. 


5 


-  VeKTvCPiL  Sherr  \*ef 
\_*TER*l  SUEUR  &,R£m 


Pitching  moment,  yawing  moment,  and  axial  longitudinal  load  are  shifted  from  the 
reference  axis  system  through  the  distances  bn,  cn  to  the  neutral  axis  prior  to 
calculation  of  the  resulting  deflections  in  the  elastic  matrix.  Similarly, vertical 
and  lateral  loads  and  rolling  moment  (fuselage  torsion)  are  shifted  from  the 
reference  axes  to  the  shear  center  axis  b#,  ce  as  shown  below: 


Therefore  F  ^  =  Fi  F*  =  F* 

F'b  s  Fz  mV  =  My  + 

Nl  3  V\«<  -  beFj+CeF'j  Myg  —  C^FX 

Similarly  the  linear  and  angular  deflections  are  transformed  to  their  respective  Set 


0<  -  C><  X'*.~byi 

fi"  ‘fi  i  *  h  •  c<*  «• 

y"  =  X  t  *  i  *  t>e 
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Xp  =  x 

^  COSC<p  +  Z  5:No<p 
Zp  =  S>»N<*.P  +  Z  COS  <XP 

FXp  -  F> 

F^p  *  F^  Cos  t^p  -4-  Fj  S.NcXp 
FeP  **p  +•  F2  cos  <xp 


*1 

cy{  p  *  ctC  i 

jS  p  =  co$  XP  -v-  'V  SlN^p 
Jf  p  =  -/5IN°<p  -*-  Y  COS  <*p 

hA  «,  p  3  NA  at 

\Afle  ~  co s  M  V  *p 

M  j  p  =  ~  Mjg  SIN  <Xp  4  VA  ^  CoS  <X  p 
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sin 


Thus,  the  triple  aatriz  product, 

1**1  fir-l  fpri m  [Elastic  I  Rotation  Translation 

L3jJ_E2jLElJ  [tutrix  J  [_HatrixJ  L  Matrix 

will  relate  the  forces  and  displacements  at  station  n  with  those  at(n  +  l)  with 
respect  to  the  principal  -  neutral  and  principal  -  elastic  axes  of  the  elastic  beam. 
Since  the  forces  and  displacements  must  be  with  respect  to  the  reference  axes  in 


order  to  continue  with  the  associated  matrix  method,  a  rotation  and  translation 


opposite  to  those  derived  previously  are  required.  Thus 
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will 


relate  the  forces  and  displacements  at  station  n  with 


those  at  station(n  +  l) 


with  respect  to  the  reference  axes. 


The  required  matrices  to  return  to  the  reference  axes  are  as  follows: 
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Hast  Matrix 


Consider  a  ma8s  M  whose  position  is  known  with  respect  to  the  reference  axes  of 
an  arbitrary  section,  n  associated  with  the  mass  are  the  inertial  properties  I  <*  , 
1 4  and  ly  about  its  centroid.  These  masses  and  inertias,  accelerated  by 
harmonic  occil lations ,  produce  inertial  loads  on  the  beam  structure. 


=  +  nnuJ,'(x^4-C^n-bVh)'-F>ltt  +  | 

]T  f  ^  -  O  =  -  F  ^  nil 

H  ‘  °  **  ^  ^  *  +b«Kh-a^n)-FlttM 

51  *  °  '  ~  c^hta  V^c.  -v  m  a^n)  b 

+  X<  -  M-<  h+( 

o  =  H^h  wwl(xh  +  c^'bVh)c-  ^cu"'(-en  +  b<*n  -  a 
+ 1/> 

5IM  Y  ~  o  ~  ^  -  ton  oo1-  ^Ob  +  c*<n  +  a-^')£i' 


Deflections  are  unchanged  across  the  mass  stations.  The  final  matrix  form  is 
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MASS  MATRIX 
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Ground  Spring  Matrix 


A  Mathematical  device  which  permits  simulation  of  ground  shake  tests,  in  which 
the  aircraft  is  suspended  from  the  roof  trusses  of  a  hangar  building,  in  the  ground 
spring  matrix.  Suppose  a  station,  n,  has  its  six  degrees  of  freedom  partially 
restrained  by  springs,  then  these  external  springs  induce  additional  loads  on  the 
vibrating  structure  in  proportion  to  the  motion  amplitude  and  spring  rates. 


f\-°*  K,(x,+cf„-bv„)- 

+  F*»+l 

°  ■  M  ^  4-K  ^  C  -v-  a  'S^c,  -  (is+bo(h-a^)  ^ 

=  O  =  (*  n  +  -  bO  ^  4  b 

~  n+l 

^  ir =  °  ~  Kjrn4  4  b  Tf*)  b  -  K^C^n- -va. 

-“v**  -HVntl 
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Concentrated  springs  bstwssn  adjacent  fuselagn  beam  sections,  such  as  the  attachment 
of  a  rotor  transmission  to  the  fuselage  structure,  are  handled  by  a  matrix.  Consider 


»>♦  i ,  Me x.  n  +i, 
n+i;-c(  o  +  i 


Since  no  force  generation  occurs  across  the  spring  joint,  the  forces  at  stations 

n  and  n  ■  1  are  equal.  Deflections  change  because  of  the  action  of  forces  on  the 

springs.  At  joint  n,  ,  s. 

I—  ix  -  O  -  rx  i\  Kx  {  Xn-  X»v»i  j 

^  Flj  -  O 

EFz-O-  Fzrv  -  -Zm-i) 

H'O'V  -  Koc 

~  K*  (p*.  -  ./Um) 

/El  My  =  O  -  M^rv.  “  Kjr  (  (fry  ~  tfm.i'i) 

Solving  these  equations  for  the  displacements  at  station  n  +  1 ,  the  concentrated 
spring  utrix  is, 
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Shift  Matrix 


To  permit  a  change  of  position  of  the  reference  axes,  a  shift  matrix  is  employed 


to  convert  loads  and  deflections  from  one  axis  orientation  to  another 


F»rv(  lv\<xn) 
x  rv^  <Xv\ 


=  K** 

CS.  p\  +  |  - 

Fy> i»* 

*  Fy> 

B  n+  l  — 

B  «. 

Fzn.l 

-  Fx. 

V  *\  1 1  s 

Yw 

tTl 

l  e  Moth  -  6*  Fx*>  +  Cofyn 

Xmi  - 

X.W  "  b  •  +  Co  On 

fUfi  n  <•* 

1  M{J* 

Co  b*  +-  Cl a  Fz  * 

Ynf-l  "  Vn  -  Co  °<^ao  Y* 

ITil'  n  r  i 

m  My* 

*  bo  F*n-  F y* 

L  n  +  >  =  Zn  +  b*  o^n  -a0  0 

Thus  the  shift  matrix 
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Bend  Matrix 


To  penile  e  change  of  direction  of  the  elastic  axis,  a  bend  matrix  is  employed 
to  convert  loads  and  deflections  from  one  axis  orientation  to  another 


In  progressing  from  left  to  right,  the  bend  angle,* is  defined  as  clockwise  from  an 
extension  of  the  initial  reference  axis  past  the  bend. 


Xx+i  ~  X*  Cos  (2. if- ft  p )  *  Sin  (zTT ~  ft p) 
Zn*i  =  Zh  COS  (ITT- pP)~  Xn.  3  l  NJ  (z  ft- ftp) 

y** i  =  i* 


<**♦1  ‘oc^css  (zir-pt>)+  (zri-p?) 

X*.4r\  =  cos  (zi7-/3p)  -  s  in  (2t r-p?) 

-  ft  rv_ 


Fxn.fi-  Pxrv  Cos  (ZTT'Pp  )+  Fzrv  ^in  (i.  TT-0  f 
Fz.*+\  -  Fzn.  COS  (Z(T-fip)~  Fx^SiN  (2ir-pp) 
fijntl  -  F'j  k 

M  «  n-n  -  M  «*v  CoS(zTF-  in  (zir-)3p) 

Mar'n*!  *  Mjfn  COS  (zV-ftp)  -  M  <*»v  5  in  (ZJh  ft p) 
Mpnvl  *  Mpiv 


Sin 


Consider  a  mass,  m,  whose  position  is  known  with  respect  to  the  axes  at  an 
arbitrary  station,  n.  Associated  with  the  mass  are  the  inertial  properties 
lot  ,  Ijt  ,  Iy  about  its  centroid.  In  ndditipn,  the  mass  is  mounted  on  six 
simple  springs  which  act  in  the  same  directions  as  the  general  displacements 
These  masses  and  inertias,  accelerated  by  harmonic  oscillations,  produce 
inertial  loads  on  the  beam  structure. 


(a)  Uncoupled  system 


H  F*,  -  °  "  Fx  n  -V  x  UJT  (  *  *  C/3  -  b  v)  - 

-  ^  ^  -  C  <*  +  a.  *)  -  F  ^  ^  ( 

5lF?  -  °  ~  F?  „  +  uo1"  (  Z  +  b  -  F*  mi 

H  M*  =  °  =  -  -  Co(+<x  CO1  (€+bcK-a^)b 

+  I  *  *  oS~  *<.  ^ 

£  =  O  =  -t-  -  b*)  C  -  VY^e  ajt(?  +  b  <^-a_p)  a 

+  V  Ufi  n+, 

^Tm^o  =  MVT,T-vw>XllUJl('>c  +  c/4-blr')b  +  irvyi  U>1(  ^  -c<*  +  a^cL 

+  1  V  H*  uji  y  ~  W 
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where  from  simple  spring-mess  theory,^/  represents  the  emplification 
factors  of  the  sys  tern 
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(b) 


The  coupled  force  relationship  (effective  mass  matrix) ,  let  mass 


properties  be  denoted  as  follows: 

m-  -  m, 

M  =  1^3 

m  fru 

I*  -  my 

Ic  -  fTU 


Let  mass  motion  be  denoted  as  follows  (with  respect  to  fuselage): 

ze  -  e, 

A  e  -  e*. 

Xe  =  e» 

y«  -  e* 

2fe  05 

®-e  fc!  b 


Let  fuselage  motion  and  external  forces  be  denoted  as  follows: 
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Let  the  mass  motion  with  respect  to  fuselage  be  represented  as  the  sum  o 


all  normal  modes 
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Where  is  the  generalized  coordinate  in  the  A  mode. 
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Substituting  in  Jyx.  equation  , 


\  W  (A->o 


ffu  f* 


-Q 


r*($ '  ’) 


•r*  / 


for  simplicity  the  equation  is  written 

00  »" 

'T'V  fj  -  -<£>a 

Since  harmonic  motion  is  assumed, 

m;  co1  +  .A  aj  f  j  uuv  -  Q  a. 

A  schematic  free  body  diagram  Gf  the  effective  mass  -  fuselage  combination  between 


where 


AMD  5l 


n-*-i 


FF-  -amu  i_*_  are 

the  six  components  of 
force  at  station  n  and 
n  +  1  respectively 


Fa  =  O  *  Fa  +-  ft  a.  -  Fa. 


or 


FJ 


r-  " 

Fa.  + 


M  i>0  ; 


t  A  rrUj  ^  u>  1 


Since  the  position  vector  of  station  n  +  1  with  respect  to  station  n  is  small,  the 
six  displacement  components  at  the  two  stations  are  equal.  Thus  the  "effective 
mass  matrix"  is  written  as  follows: 
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COUPLED  SPRUNG  MASS  MATRIX 


Force  Matrix 


A  mathematical  device  which  permits  simulation  of  a  forcing  function  placed 
anywhere  along  the  fuselage  of  an  aircraft  is  the  force  matrix.  Suppose  a 
station,  n,  is  subjected  to  a  forcing  function  having  six  components,  then  these 
external  forces  induce  additional  loads  on  the  vibrating  structure. 


LFx-O'  F*„-  F*  -  Fw, 

^  -  O  - 

F  ?r  =  O  ~  + 

ZZMo^-o  =  M**  -  C?^  -V-  b?z-  tAKn+l  +  M* 

I>,=  o  =  H^+cF^  V\p 

MT  -  O  =  My  -  bfx  -v-Q-fvj  —  Myn+,  -*r  w\y 
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is  of  the  form 


In  general,  the  forcing  function,  Fij 

r-Fo  +  F-,  COSY„  +- 

where  vf»n  is  harmonic  azimuth  angle 

Sine.  the  forcing  function  .u.t  be  edded  to  the  lnertl.1  nnd  eln.tlc  force. 

■PPlied  to  the  fue.l.g.,  the  force,  ere  b.ndled  by  eddlng  .  Uth  to.  end  colon 
to  the  previous  12  x  12  matrix. 

This  13  x  13  force  matrix  requires  that  all  the  previously  derived  matrices 
be  padded  from  12  x  12  matrices  to  13  x  13. 
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FORCE  MATRIX 
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Unit  Forces  (lb.)  F&*  ,  Fc,^  ,  Sin  Force  Coefficients  (#) 

Niu(5  Unit  Moaents  (In.#)  Ms<x,  Sin  Moment  Coefficients  (In.#) 

'  c  Cos  Force  Coefficients  (#)  >Y„  Haraonic  Azimuth  (deg.) 

m’  m  w  Cos  Moment  Coefficients  (In.#)  cl.  b«  CT  Lo°g>  lat.  &  Vert.  Force  Locations 


General  Matrix 


The  general  matrix  shown  below  was  placed  in  the  computer  program  to  permit  the 
formation  of  any  additional  matrices  needed  in  the  coupled  vertical -lateral 
matrix  analysis. 
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Having  the  racking  coefficients  for  relative  translations  in  the  vertical  direction, 
enables  calculation  of  the  relative  vertical  translation  at  any  point  in  the  fuselage 
and  the  relative  rotation  about  the  longitudinal  axis 


w 


Zn+i  -~2_n  =  04  Moc-  n  +  (°c "+• 1  _0<:  n)  u/| 


Replacing  the  relative  rotation  in  the  above  equation 

Z  nti  -  Zn  =  ('  ^  )  0CCX.  Mcc.n  *  ^ 

Likewise,  using  the  lateral  racking  coefficients 

0  b  =<  M  a.  n 


<pd  °C  M  «n 


Placing  the  relative  rotation  in  the  above  equating 


V,.,  -  Yn  *  (l-  (p  b<K  M  <\r)  +  _H?-  (p  d  x  M pt n 
\  TT  )  H 
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Thus,  if  ths  squat ions  for  vertical  and  lateral  translation  and  longitudinal 
rotation  are  added  to  the  equations  comprising  the  existing  elastic  matrix,  the 
result  is  an  elastic  matrix  which  includes  frame  racking. 
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The  frame  racking  matrix  in  the  lateral  and  vertical  direction  is  used  with 
neutral  axis  and  shear  center  axis  shift  matrices  developed  earlier.  The  sub¬ 
divided  elastic  matrix  can  be  written  as. 


r  ~ 

r 

—  — 

r  1 

r 

H 

1 _ 

1 

1 _ 

1 

1 _ 

Input  required  for  E^,  E2,  E3  matrices  shown  on  the  following  pages. 

bQ,  Dist.  from  the  ref.  axis  to  the  neutral  axis  in  the  y  direction  (in.) 

Cn,  Dist.  from  the  ref.  axis  to  the  neutral  axis  in  the  z  direction  (in.) 

2 

Ax,  Axial  Compression  area  (in  ) 

9 

Ay,  effective  shear  area  in  the  y  direction  (in  ) 

2 

Ax,  effective  shear  area  in  the  z  direction  (in  ) 

4 

1^,  Bending  stiffness  about  the  y  axis  (in  ) 

I£,  Bending  stiffness  about  the  z  axis  (in^) 

2 

E,  Bending  modulus  (#/in  ) 

2 

G,  Shear  modulus  (#/in  ) 

L,  Elastic  bay  length  (In) 

fia.*.,  Linear  deflection  at  the  right  side  of  the  frame  from  a  unit  moment  (1/#) 
Linear  deflection  at  the  top  section  of  the  frame  from  a  unit  moment  ( J  / #) 
(&:ocj  Linear  deflection  at  the  left  side  of  the  frame  from  a  unit  moment  (1/#) 

(fid  06 ,  Linear  deflection  at  the  bottom  section  of  the  frame  from  a  unit  moment  (1/#) 
W,  Lateral  distance  between  longerons  (in.) 

W^,  Lateral  distance  from  the  reference  axis  to  the  right  longeron  (in.) 

H,  Vertical  distance  between  longerons  (in.) 

H2,  Vertical  distance  from  the  reference  axis  to  the  upper  longeron  (in.) 
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effective  wing  matrix 

Inclusion  of  the  floating  fuel  wings  in  the  representation  requires  an 
effective  wing  matrix  which  transmit  the  inertial  loads  of  the  wing  to 
the  fuselage  beam  station.  Consider  the  plan  view  of  the  fuselage-wing 
as  shown  below. 


|  FWD 


*yi  >  ^  1  j  <X ^  |5,  ^  GENERAL  DISPLACEMENTS  OF  FUSELAGE 
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^  +  Q-t  tf 

2:  c  —  2.  +  i  c<  —  aL  |3  -h  2?  ^  — 

Summing  along  the  wing,  the  kinetic  energy  of  the  right  wing  is, 


T  *  T  Im; 


+  +  o.L  y ) 


t  +  Cz  4  ^  +  ZL  -  o-i. 

+iZ  {  l,c-5 * «S  *  v,  ( p +  po1*-  c  <  -  -sJ  ) 
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Expressing  the  right  wing  generalized  loads  in  matrix  form  as  a  function 
of  displacement 


co^Ero  i. 

co  E-r^a^ 

r  otuS* 

r-<7eCt  tu* 

“C01£  rr)La.L 
rcrto,co* 

kfcS-0 

k(jg-0 

to*Erpi>4ri 

ufS'miir'+O 

.V 

nci-') 

k(*H 

hm-t) 

*.r-  2 

rC>ara  to1 

k(3i-0 

y  or/- 

kcS-0 

m 

-  ojrJ7  m  ;  Jbi 

In  terms  of  the  effective  spring,  Ke  in  each  mode,  the  potential 
energy  of  the  right  wing  is 

tt») 

Ke  Hs 

S-l 

where,  (*) 

Ke=  Is  "s 

Applying  Lagrange ' 8  equation,  the  right  wing  equations  of  motion  are, 

=  F * 

©  o 

1  (Em;.)  2  +  +  (-Eitu  0-0/3° 

-*•  H,  »  F, 
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+  (r£nuajiri)  p 

+  [EwA(Z?'a'<.^)  +  II4(i^<i.]  Hs  =  M<* 

/3  [Owia-c‘+  ifO]p  +(-z»«tO'i.)z 

+  [”Z  m.a,,  ft  )+  Hs  = 

tf  [XXmia*  +miAt  +  Iv„)j  +(£nnal)<^  =  M* 

Hs  V.«? *-*■  H°s  +  KsHs 

+  [I  mi  (**-«*  rf*}] 


Similarly,  the  left  wing  generalized  loads  are  obtained  and  then,  combined 
with  the  right  wing  to  obtain  the  following  expressions. 


Maw 


-  Maw 


Combining  the  deflections  with  the  above  load  expression,  the  effective 
wing  matrix  is. 
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EFFECTIVE  WING  MATRIX 
(H-21  RANGE  EXTENSION) 


APPENDIX  B 

METHOD  OF  SOLUTION  FOR  COUPLED  MATRIX  PROGRAM 
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METHOD  OF  SOLUTION  FOR  COUPLED  MATRIX  PROGRAM 


* 


Method  of  Solution  -  Free  Vibration 

These  matrices  ere  now  "building  blocks"  from  which  e  dynamic  representation 
of  the  helicopter  fuselage  may  be  constructed. 

The  boundary  conditions  for  such  a  representative  helicopter  fuselage  are  that 
the  loads  Fx,  Fy,  Fz,  M^,  and  My  just  before  the  forward  rotor  mess  and 
Just  after  the  aft  rotor  mass  are  zero,  that  is  a  "free-free"  beam.  Even 
ground  springs  do  not  affect  the  boundary  conditions,  because  they  are  taken 
to  act  on  the  fuselage  side  before  the  boundaries  are  reached.  Matrix  multi¬ 
plications  from  the  forward  to  aft  rotor  are  conducted  repeatedly  for  trial 
frequencies,  until  the  boundary  conditions  are  met.  Note  that  since  the  for¬ 
ward  end  boundary  conditions,  Fx  ■  Fy  -  Fx  ■  Mot  “  ■  Mjj  -  0,  appear  in  the 

single  column  end  matrix,  and  since  only  these  forces  need  be  determined  for 
comparison  with  the  aft  end  boundary  conditions,  it  is  then  necessary  to  use 
only  the  terms  shown  on  the  next  page  in  the  final  collapsed  matrix  for  the 
frequency  calculation. 
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Repeated  trials  of  frequency,  w,  are  raade,  and  each  one  of  the  thirty-six 
terms  shown  above  are  obtained  numerically.  To  meet  the  boundary  requirement, 
the  residual,  A,  of  the  matrix  array  must  be  zero. 


a,5 

a... 
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nt 
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In  practice,  the  residual  may  be  plotted  against  the  frequency  trials,  and  when 


A 


a  aero  intersection  occurs,  natural  frequencies  uj,  ,  uJ„a,'e  determined- 
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For  each  natural  mode,  It  la  desirable  to  normalize  In  terms  of  the  largest 
relative  linear  deflection  value  which  exist  at  the  ends  of  the  beam.  Now, 
calculating  the  relative  linear  deflections  at  both  ends  of  the  beam. 

Rewriting  the  collapsed  matrix  with  the  zero  columns  removed. 


Aft  ib  5  *  .  >  i‘ 

Using  th'!  aft  rotor  load  equations, 


V  >J  s/t»r 
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FORCED  VIBRATION 


Forced  mode  shapes  can  be  obtained  by  collapsing  the  matrix  array  representing  the 
fuselage  into  a  single  matrix  considering  the  unknown  frequency, co  as  the  frequency 
of  forced  excitation.  Including  the  known  magnitude  of  the  shaking  loads  in  the 
fuselage  system  the  single  matrix  for  the  entfrf  fuselage  becomes. 


Aft  Rotor  Fwc1  Rotor 


E 


1) 

To  solve  for  tho  unknown  quantities  at  the  forward  rotor,  consider  the  boundary 
conditions  of  the  forward  and  aft  rotor  as  part  of  its  corresponding  coluun  matrix 
and  then,  simplify  the  solution.  As  an  example,  consider  a  free-free  beam;  the 
boundary  conditions  are  defined  by  making  all  external  loads  equal  to  zero.  There¬ 
fore,  the  aMtrix  can  be  written  as. 
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How,  the  general  solution  considers  only  the  equation  which  are  equal  to  zero.  There¬ 
fore  the  boundary  conditions  can  be  represented  by  six  linear  equations.  The  matrix 
for  the  free-free  beam  becomes, 
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Th«  6x6  boundary  condition  matrix  may  be  obtained  directly  by  considering 
only  the  elements  common  to  rows  of  zero  boundary  at  the  aft  rotor  and  to 
columns  of  non-zero  boundary  condition  at  fwd  rotor.  Considering  the  use  of 
this  procedure  in  the  free-free  beam  cast 


Continuing,  the  matrix  con  bo  rovrltton  In  torso  of  Q  which  io  o  known 
toIuo  for  tho  fsno-froo  boas. 
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Tho  oquotlono  con  now  bo  solved  for  tho  unknown  voluoo  by  inversion  of 
tho  square  matrix.  Tho  solution  for  tho  freo-frea  boos  con  be 
written  os, 
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Following  avoluotlon  of  the  fwd  rotor  boundary  conditions  tho  forced 

sodo  shape  at  aoeh  station  io  obtained  by  progressive  matrix 

suit ip 11 cot ion  and  print-out  storting  with  tho  known  colusn  elements. 
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